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1 Introduction 

B. H. Neumann [N] introduced outer billiards in the late 1950s and J. Moser 
[Ml] popularized the system in the 1970s as a toy model for celestial mechan- 
ics. Outer billiards is a discrete self-map of R 2 - P, where P is a bounded 
convex planar set as in Figure 1.1 below. Given pi e R 2 — P, one defines 
p 2 so that the segment pvpi is tangent to P at its midpoint and P lies to 
the right of the ray p\p\- The map p\ — > p 2 is called the outer billiards map. 
The map is almost everywhere defined and invertible. See [Tl] for a survey 
of outer billiards. 





Figure 1.1: Outer billiards relative to P. 
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Usually, the main object of study is the orbit {p n } of the point p\. These 
orbits are both complicated and interesting. When P is a convex polygon, 
the orbits frequently have a fractal structure. Aside from a few cases, the 
structure of the orbits is not yet well understood. 

It turns out that it is sometimes productive to study not the orbit it- 
self but rather a certain acceleration of the orbit. Roughly speaking, the 
acceleration we have in mind amounts to considering the first return map to 
a certain infinite strip in the plane. We call this acceleration the pinwheel 
dynamics, because it is based on the pinwheel map that we studied in [S3]. 
We give a full account of the pinwheel map in §2.21 and recall some of the 
results of [S3] in §2.51 We point out, however, that this paper does not rely 
on any of the results in [S3]. The work here is self-contained. 

The first purpose of this paper is to introduce an object called the arith- 
metic graph. The arithmetic graph is a polygonal path T(P, p) C R n that 
one associates to the pinwheel dynamics of p = p\. One can view this graph 
as a geometric incarnation of the symbolic coding of the pinwheel dynamics. 
The arithmetic graph is quite similar to the lattice paths studied by Vivaldi 
et. al. [V] in connection with interval exchange transformations. We will 
define the arithmetic graph in §2.3[ right after defining the pinwheel map. 

In case P is a kite - i.e. a convex quadrilateral having a diagonal that is 
a line of symmetry - we have studied the arithmetic graph in great detail. 
In [SI] we analyzed the graph T(P,p), where P is the Penrose kite and p is a 
specially chosen point. For this pair, we showed that T(P,p) is coarsely self- 
similar, in the sense that a certain rescaled limit of T(P,p) is a self-similar 
fractal curve. See §4. II for a formal definition of what we mean by a rescaled 
limit. 

We used the coarse self-similarity of the arithmetic graph in this case to 
conclude that the orbit of p is unbounded relative to P. This provided the 
first example of an outer billiards system with unbounded orbits. In [S2] we 
studied the arithmetic graphs relative to kites in general, and showed that 
outer billiards has unbounded orbits when defined relative to any irrational 
kite. In §3.11 we will show some of the nicest pictures of arithmetic graphs 
associated to kites. 

The irrational kites are the only known examples of polygons for which 
outer billiards has unbounded orbits. Since the arithmetic graph turns out 

1 It is worth mentioning that Dolgopyat and Fayad show in [DF] that outer billiards 
also has unbounded orbits when defined relative to a semi-disk. 
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to be a decisive tool for establishing the unboundedness results for kites, we 
think that it will also be a useful tool for polygonal outer billiards more gen- 
erally. Compare Theorem 12.31 (a result quoted from [S3]) and the discussion 
following it. 

Even in the case when all the orbits are known to be bounded, as they are 
for so-called quasi-rational polygons - see [VS], [K], [GS] - we think that 
the arithmetic graph should shed light on the dynamics. (See the end of §2.51 
for a definition of quasi-rational.) In particular, the arithmetic graph sheds 
light on the dynamics relative to the regular polygons. The regular polygons 
are nice examples of quasi-rational polygons. 

The vertices of a regular polygon have coordinates that lie in a cyclotomic 
field. When p lies in this same field, one can view certain projections of V(p) 
as "Galois conjugates" of the orbits. In somewhat the same way that the full 
list of Galois conjugates of an algebraic number sheds light on the number, 
the arithmetic graph sheds light on the orbits associated to a polygon with 
algebraic vertices. 

This brings us to the main purpose of our paper. We will concentrate on 
outer billiards for the regular octagon, a polygon that we scale so its vertices 
are the 8th roots of unity. The outer billiards dynamics for the regular 
octagon are well understood from certain points of view, but we will show 
that the study of the arithmetic graphs in the regular octagon case leads to 
a big surprise: There is a fractal loop T in R 4 which has the following three 
properties. 

• One planar projection 7r 3 of T is an embedded fractal loop reminiscent 
of the Koch snowflake. See Figure 1.2 below. We call this fractal the 
snow flake. We give a precise definition in §4.21 

• Another planar projection tt 2 of T is a fractal set reminiscent of the 
Sierpinski carpet. See Figure 1.3 below for an approximate picture. 
We call this fractal the carpet. We give a precise definition in §4.41 

• T is the rescaled limit of the arithmetic graphs associated to any se- 
quence of odd periodic orbits. We give a precise definition of rescaled 
limit in §4.11 

It turns out that all the periodic orbits intersect a certain strip in 3 fe 
points for some integer k. See §7.41 We call the orbit odd if k is odd. In 
Figure 1.2, the odd orbits are arithmetically closed orbits are lightly colored. 
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What we actually prove is slightly weaker than what we have just said. For 
technical reasons, we only consider odd orbits that lie outside the first layer 
of big grey octagons in Figure 1.3. 




Figure 1.2: The snowflake and the carpet 




Figure 1.3: Periodic orbits associated to the octagon 

By considering the correspondence between our two fractals, we produce 
a surjective continuous map from the snowflake to the carpet. This map is 
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"symmetry respecting" in a certain formal sense that we discuss in §H See 
§4.61 for a description of the map from the snowflake to the carpet. What 
is surprising is that this canonical map between two seemingly unrelated 
fractals actually arises "naturally" , in connection with outer billiards on the 
octagon. 

For the sake of comparison, we now discuss the situation for other regular 
polygons. For n = 3, 4, 6, outer billiards on the regular n-gon is rather trivial 
and easy to describe. In all these cases, there is a familiar periodic tiling of 
the plane, such that the dynamics permutes the tiles. In the case n = 4 we 
get the usual square tiling. In the cases n = 3,6 we get a tiling by hexagons 
and triangles. 

For the case n = 5, Tabachnikov [T2] analyzes the situation in detail. 
The recent paper [BC] builds on the work in [T2] and describes the sym- 
bolic dynamics for the periodic orbits in great detail. We are interested in 
something different, namely the pinwheel dynamics, but the two are closely 
related. In the case of the regular pentagon, our construction produces a pair 
of embedded fractal curves, both akin to the Koch snowflake, that arise as 
geometric limits of arithmetic graphs associated to periodic orbits. See §3.21 
For the sake of brevity, we do not present proofs of the statements we make 
for the regular pentagon. We are mainly interested in the octagon, and the 
analysis we give in the case of the octagon could be fairly easily redone for 
the pentagon. 

The four cases n = 5, 10, 8, 12 are all quite similar to each other, and 
possible to understand completely. In these cases, outer billiards on the 
regular n-gon has an efficient renormalization scheme that allows one to 
give (at least in principle) a complete description of what is going on. The 
special properties of these values of n is that the corresponding nth roots of 
unity lie in quadratic number fields. In the case n = 10, we get the same 
two snowflakes as in the case n = 5. We leave the case n = 12, which 
is mildly more complicated than any of the other quadratic cases, to the 
experimentally minded reader. 

For any other positive integer n, outer billiards on the regular n-gon is 
not well understood at all. For the sake of comparison, we will draw some 
pictures of arithmetic graphs associated to the regular 7-gon. See §3.41 In 
this case, the pictures reveal an explosion of complexity in a geometrically 
striking way. Our 4 figures in §3.41 give only the faintest hint of the complex- 
ity. 
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Here is the plan of the paper. 

• Basic Definitions: In §2 we will define the arithmetic graph, and give 
explicit formulas in the case of the regular octagon. 

• A Picture Gallery: In §3 we draw the arithmetic graphs for many 
examples, including the regular octagon. After reading §3.31 the reader 
should have a good intuitive idea of what our main theorem will say, 
even though we defer the formal statement until §5.11 

• Carpet and Snowflake: In §4 we define the snowflake and carpet 
fractals precisely, and discuss the relationship between them. 

• Reduction to a Small Region: We state our main result, Theorem 
I5.2[ at the beginning of §5. In §5 we reduce Theorem 15.21 to the study 
of periodic points in a very small region Ri. 

• Toy Example of Renormalization: In §6 we study a simple and 
well-known polygon exchange map that arises in connection with outer 
billiards on the regular octagon. The material in this chapter is not 
strictly necessary for our formal proof, but the system we study here is 
closely related to the one that comes from the pinwheel map. 

• Pinwheel Dynamics: In §7, we describe the dynamics of the pinwheel 
map on the region R\. Actually, we reduce an even smaller region, R, 
which is the top half of R\. The dynamical system on R is the main 
dynamical system of interest to us. Both the pinwheel map and the 
system from §6 have a period 3 renormalization that is the key to their 
analysis. 

• Substitution Scheme: In §8 we study the arithmetic graphs pro- 
duced by the dynamics on R. Using the renormalization scheme, we 
describe an efficient 2-part method for generating the projections of the 
arithmetic graph of interest to us. The first part of the method is com- 
binatorial, and involves a substitution scheme for numerical sequences. 
The second part is geometrical, where we replace each number in the 
created string by a certain vector. We call this second half the vector 
assignment. 
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• Fixed Point of Renormalization: In §9 we modify the construction 
presented in §8, keeping the combinatorial part the same but changing 
the vector assignment. We see that the substitution scheme produces a 
kind of renormalization operator defined on the vector space of all pos- 
sible vector assignments. This fixed point corresponds to the Perron- 
Frobenius eigenvector. By taking the fixed point of our operator as the 
new vector assignment, we produce a nicer family of curves that has 
the same rescaled limit. 

• Self-Similar Pattern Matching: In §10 we check by a combina- 
tion of direct calculation and induction that our improved curves from 
§9 have the snowflake and the carpet as their limit. The key idea is 
that the improved curves and the fractals from §4 are self-similar in 
compatible ways. The final analysis completes the proof of the Main 
Result. 

The reader might be interested to know that we have made two java pro- 
grams, OctoMap 1 and OctoMap 2, which illustrate the mathematics in this 
paper. OctoMap 1 illustrates the snowflake and carpet fractals in great de- 
tail, and gives an interactive demonstration of the material in §4. OctoMap 
2, a much more extensive program, starts with the dynamical system pro- 
duced in §7 and illustrates all the key ideas that go into the proof of the 
Main Theorem. In particular, one can use OctoMap 2 to survey all the com- 
putations we describe in §7-10. One can find our applets, respectively, at the 
following address: 

• http://www.math.brown.edu/~res/Java/OctoMap/Main.html 

• http:/ /www. math. brown.edu/~res/Java/OctoMap2/Main. html 

We strongly encourage the reader to look at these applets while reading this 
paper. They relate to this paper the same way that a cooked meal relates to 
a recipe. 

I'd like to thank Gordon Hughes and Sergei Tabachnikov about interesting 
and inspiring conversations about outer billiards on regular polygons. 
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2 The Arithmetic Graph 



2.1 Strip Maps 

Let E C R 2 be an infinite strip. Let V be a vector whose tail lies on one 
edge of E and whose head lies on the other. All this is shown in Figure 2.1. 




Figure 2.1: A strip map 



To (E, V) we associate a strip map f : R 2 — >■ E, defined by /(p) = p+nV, 
where n E Z is chosen so that p + nV G E. In the example shown, we have 
n = 2. 

Remarks: 

(i) / is not everywhere defined. It is not defined on a certain countable family 
of lines that are parallel to E. In particular / is not defined on the boundary 
of E. 

(ii) The pair (E, — V) gives rise to the same map as the pair (E, V). 

(iii) In case E is the horizontain strip bounded by the lines y — and y — 1, 
and V = (0, 1), the map map / has the formula f(x,y) = (x, [y]), where [y] 
is the fractional part of y. In this case, / is not defined on the horizontal 
lines of integer height. Any other strip map is conjugate to this one by some 
affine transformation. 
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2.2 The Pinwheel Map 

Let P be a convex n-gon. To P we associate n special pinwheel strip. Each 
pinwheel strip E is such that one component of <9E contains an edge e of P, 
and vertices of P farthest from this boundary component lie on the centerline 
of E. See Figure 2.2 




Figure 2.2: A pinwheel strip associated to a regular octagon 

We orient the boundaries of a pinwheel strip E so that a person walking 
along a boundary component would see E on the left. Thus, the two boundary 
components are oriented in opposite directions. Say that a pointed strip is a 
strip, together with a choice of boundary component. We think of a pointed 
strip as pointing in the direction of the orientation on its preferred boundary 
component. To each n-gon we associate In pointed strips, 2 per pinwheel 
strip. We denote pointed strips by a pair (E,L), where L is a boundary 
component of the strip E. 

To the pair (E, L) we associate a vector V as follows. We choose e > 
and let p be a point which is e units from L, outside of E, and 1/e units 
away from the origin. Of the two possible locations for p that the above 
conditions determine, we choose the one toward which L points. Figure 2.2 
shows what we mean. We then let V = 4> 2 {p) — p. Here 2 is the square of 
the outer billiards map. Our definition is independent of e, provided that e 
is sufficiently small. 
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In general, V = 2(w — v), where v is the vertex bisecting the segment 
joining p to 4>(p) and w is the vertex bisecting the segment joining <p(p) to 
4> 2 {p). It is not hard to check that £ and V are related as in §2.11 




12 13 

Figure 2.3: The pointed strips associated to a regular octagon 

Thus, we have associated to P a total of 2n triples (S fc ,L fc ,Vfc). Here 
(£fc,Lfc) is a pointed pinwheel strip and Vk is the associated vector. We 
cyclically order these In triples in the following way. We choose a large 
circle C centered at the origin, and orient C counterclockwise. Each pair 
(Sfc, Lfe) determines a point Cfc G C, namely the intersection point of PI P 
towards which points. We order our triples so that the points ci, ...,C2„ 
occur in counterclockwise order along C. Figure 2.3 shows these points, in 
case P is the regular octagon. In Figure 2.3 we highlight the pointed strips 
E 15 and S = E 16 . 
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Let fk be the strip map associated to the pair (S fe , Vk). We say that the 
pinwheel map is the composition 

$ : I'2n o . . . o f : S S (1) 

We say that p G S is preferred if, for all x E P, the vector p — a; has pos- 
itive dot product with vectors parallel to Lq. We call the subset of preferred 
points the preferred half of So and denote it by Hq. We usually restrict <£> to 
H Q . Note that <£> need not carry Hq to Hq. However, if p G Hq is sufficiently 
far from S - H then $(p) G i? - 

Remark: For our purposes here, the boundary of Hq (that separates Hq 
from So — Hq) is not important. What is important is the end of So that 
Hq determines. Below, and just for the case of the regular octagon, we will 
replace H by a similar set, namely Sq below, that is invariant under $. 

Let <j) denote the outer billiards map and let <fi 2 be the square of the outer 
billiards map. Let O+(0 2 ,p) denote the forwards 2 -orbit of p. 

Lemma 2.1 The following is true for all points p G Hq that lie sufficiently 
far from P. Let q be the first point in O + (0 2 ,p) that lies in Hq. Then 
g = $(p). 

Proof: For each point p G R 2 — P for which 2 is well-defined, there is a 
vector V(p) such that <p 2 (p) —p = V(p). The function p — > V(p) is locally 
constant. Let A be the disk bounded by the large circle C, as in Figure 2.2. 
Let Wk be the open acute wedge shaped region of R 2 — A bounded by the 
lines L k -i and L k . The vector V k has the property that V(p) = Vk for at 
least some points p G W^. But one can check easily that 2 is defined on all 
of Wk provided that the circle C is taken large enough. Hence V(p) = Vk for 
all p G W k - 

Starting at p = p G W\ fl H , the successive points in O + (0 2 ,p o ) have 
the form p + mV\ for m = 1, 2, 3.... This continues until we reach a point 
Pi — Po + ™>iVi ^ ^l- But then p 1 = fi(po). Starting at p\ G W 2 H S 1; the 
successive points in + (<j) 2 ,pi) have the form pi + mVi for m = 1,2,3.... This 
continues until we reach a point p 2 = pi + m-iVi G S 2 . But then p 2 = /^(pi)- 
Continuing in this way, we get p n G S . But p n and pq lie in opposite com- 
ponents of S — A. We have gone halfway around. Continuing the process, 
we finally arrive at p 2n = q = ^(p ) E H , as claimed. 6 
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2.3 The Arithmetic Graph 



Let p = po G So be a point. As in the proof of Lemma |2~TI there are integers 
mi, rri2n such that 

Pk ■= fk{Pk-x) = Pk-i + m k V k e Sjfc. (2) 

Let vi,...,v n be the vertices of P, ordered clockwise around P. Let 
ei,...,e n be the standard basis of R n . For each vector V — Vu there are 
indices i and j such that V k = 2(t>j — Vj). Of course, i and j depend on k, 
but we are suppressing this from our notation. We define 

V = 2{e l -e j )- fc = l,...,2n. (3) 

We define 

2n 

$(p) = ^m fc i4; $ fc+1 (p) = i(i fe (p)). (4) 

fc=i 

We define T(P,p) to be the polygonal path in R n that starts at and has 
consecutive vertices $(p), $ 2 (p), etc. 

Remarks: 

(i) It might seem at first that $ is somehow a "lift" of $ to R n . This is only 
partially true. What is true is that $ fc is a map from So to R n . This map is 
only defined on points where $ fc is well defined. 

(ii) By construction the arithmetic graph is a lattice path in that its vertices 
lie in Z n . As remarked in the introduction, one should compare the lattice 
paths of [V]. 

(iii) There is a constant C, depending only on P, such that \mk+ n — mk\ < C'. 
The reason is that the successive points pi, ...,P2n are the vertices of a convex 
polygon that is centrally symmetric to within a bounded error. Hence, the 
edges of T(P,p) have length at most C, for some constant C depending only 
on P. This property is important when we take rescaled limits of the graph. 

(iv) There is an obvious projection 7r : R n — > R 2 given by n(ek) = v k - By 
construction, p + n(T(P,p)), meaning the translation of ir(T(P,p)) by p, is 
exactly the forwards $ orbit of p. What makes T(P,p) interesting is that 
sometimes other projections are much more revealing. We will illustrate this 
in the next chapter with many examples. 
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2.4 Formulas for the Octagon 

Since we are concentrating on the regular octagon, it seems worthwhile giving 
explicit formulas in this case. Recall that the vertices of the regular octagon 
are u k for k — 0, 7. We let (ij) denote the line through u l and u k , oriented 
from % to j. We let k(ij) denote the image of (ij) after we apply reflection 
in u k . Finally, we let [ij] denote the vector that points from uu l to 2a;- 7 ' — a^. 
The vector —[ij] is literally the negative of [ij]. 




Figure 2.4: The triple (E , L , V ). 



Definition: The shaded region in Figure 2.4 is the subset of S to the right 
of the line y = 2 + \/2. We let denote this set. 

To illustrate our notation by way of example, The triple (E ,L ,Vo) is 
specified by 

: 6(23) (23) - [26]. (5) 

We first list L , then the other component of E , and then V . See Figure 
2.4. 
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The listing for (E 2 , L 2 , V 2 ) is obtained from the one for (E , L Q , V ) simply 
by incrementing all the indices by 1. That is 



2 : 



7(34) (34) 



[37]. 



And so on. 




Figure 2.5: The triple (£ 1; L x , V~i). 
The triple (Si, Li, V\) is specified by 

1 : (67) 3(67) [63] 



(6) 



See Figure 3.4. 

The listing for (E3, L3, V3) is obtained from the one for (Ei, L\, V\) simply 
by incrementing all the indices by 1. That is 



3 : 



(70) 4(70) [74] 



And so on. 
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2.5 Discussion 



Lemma 12.11 gives a precise relationship between the pinwheel map and the 
outer billiards map for points in Hq that are far from P. The relationship 
between <fr 2 and $ for points near P is rather subtle, and it was the purpose 
of the paper [S3] to explain the relationship. Here we summarize some of 
the results from [S3]. These results are not needed for this paper, but they 
are worth knowing. 

Let 0+($,p) denote the forward $-orbit of p e Hq. Here is a consequence 
of our work in [S3] 

Theorem 2.2 Suppose that all our constructions are made with respect to 
a convex polygon with no parallel sides. Then, the following holds for any 
sufficiently large disk A centered at the origin. Let p G X — A. Then 
+ ((f> 2 ,p) returns to H — A if and only if + ($,p) returns to H — A, and 
the two points of return are the same. 

There are two subtle points to Theorem 12.21 First, the forward orbits 
of p, in either case, might wind many times around P before returning to 
H — A . Second, Lemma I2TT1 only applies to points of H that are sufficiently 
far from P. After our orbits exit H — A, then might wander quite close to 
P and only "come back out" much later on. However, magically, they "come 
back out" in exactly the same way. One consequence of Theorem 12.21 is 

Theorem 2.3 Suppose that all our constructions are made with respect to 
a convex polygon with no parallel sides. The pinwheel map $ has unbounded 
orbits relative to P if and only if outer billiards has unbounded orbits relative 
to P. 

We only proved these results for polygons without parallel sides because 
we wanted to avoid certain technical complications. We fully expect that 
the same results hold for all convex polygons. However, we have not yet 
worked out the details. For the present paper, which deals (rigorously) just 
with the regular octagon, we will prove something stronger than Theorem 
12.21 Namely, in §5.31 we prove 

Lemma 2.4 (Invariance) For k = 1,2,3... let £§ C So be the subset con- 
sisting of points lying to the right of the line x = k(l + \/2) . Then 
^-invariant for all k > 1 
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The argument in Lemma 12.11 works directly for all orbits in Ej. Ac- 
cordindly, for points in Ej, we can be sure that the pinwheel dynamics and 
the first return dynamics of the square outer billiards map coincide. 

The Invariance Lemma implies that all forward orbits of both <p and <3? are 
bounded. By symmetry, the same goes for the backwards orbits. Hence, the 
Invariance Lemma implies Theorem 12.31 for the regular octagon. Compare 

HI 

The Invariance Lemma is a special case of a general result concerning 
quasi-rational polygons. The polygon P is called quasi-rational if it may be 
scaled so that 

area(E fc n E fc+1 ) G Z U oo 

for all k = l,...,2n. Here indices are taken mod 2n, as usual. In case no 
sides of P are parallel, the above areas are all finite. For regular polygons, 
one can scale so that all the finite areas are 1. Hence, regular polygons are 
quasi-rational. Likewise, polygons with rational vertices are quasi-rational. 

As we mentioned in the introduction, it is proved in [VS], [K], and [GS] 
that all outer billiards orbits are bounded for a quasi-rational polygon. In 
[S3] we give a self-contained proof of this result, in case P has no parallel 
sides. 
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3 Examples of Arithmetic Graphs 



3.1 Rational Kites 

As in [S2] we let K(A) be the kite with vertices 

(-1,0); (0,1); (0,-1); (AO); A e (0, 1). (7) 

The case A = 1 corresponds to the square, a trivial case we ignore. K(A) is 
called (ir)rational iff A is (ir)rational. 

When A = p/q, we call the point (1/5,1) the fundamental point, for 
reasons we explain in great detail in [S2] and briefly as follows. The orbit 
of any point (x, 1) has the same combinatorics as the point (1/q, 1) for any 
< x < 2/q. For this reason, the point (l/q,l) is a representative of the 
dynamics of points "arbitrarily close" to the top vertex of K(p/q) and on the 
same horizontal line as the vertex. 

We will consider the arithmetic graphs 

T p/q = T(K(p/q),(l/q,l)) (8) 

for various choices of p/q. We will draw a certain projection of this graph into 
the plane. We will not specify the projection we use explicitly because it is 
tedious to do so. The fact is that T p / q lies in a thin neighborhood of a 2-plane 
in R 4 and a random projection will produce pictures very similar to ours. 
Indeed, any projection looks about the same, up to an afhne transformation 
ofR 2 . 

Consider the rational sequence {a n } which starts 0, 1 and obeys the rule 
a n = 4a n _i + a n _2- This sequence starts out 

0,1,4,17, 72,305,... 

We have 

lim a n -i/a n = \fh — 2 = <p~ 3 . 

n^oo 

Here is the golden ratio. The kite K((j)~ 3 ) is afnnely equivalent to the 
Penrose kite. The first four quotients are 

1/4 4/17 17/72 72/305... 

Figure 3.1 shows Ta for the first four of these quotients. 
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Figure 3.1: T A for A = 1/4, 4/17, 17/72, 72/305. 



We want to emphasize that our projection is such that the vertices of 
these polygons all lie in Z 2 . Were the polygons drawn to scale, each one 
would be about 3 4.23 times as large as the predecessor. However, when 
we rescale them so that they are about the same size, a fractal structure 
emerges. Later in the paper, we will formalize what we mean to taking the 
rescaled limit of a sequence of arithmetic graphs, but we hope that the above 
informal discussion makes the general idea clear. 

Another nice sequence of pictures is given by the quotients of the rational 
sequence {a n } where ao = 1 and a\ — 2 and a n = 2a n _i + a n _ 2 . In this case 
lima n _i/a n = \/2 — 1. One of the quotients is 169/408. Figure 3.2 shows 
the picture of the corresponding arithmetic graph. One can see the fractal 
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structure emerging just from this single polygon. 




Figure 3.2: T A for A = 169/408. 



There is one interesting feature of K p / q that one sees almost immediately. 
If pq is even, then K p / q is a closed embedded loop. If pq is odd, then K p / q 
is an open embedded polygonal curve. The difference derives from the fact 
that the orbit of (1/q, 1) is stable when pq is even and unstable when pq is 
odd. By this we mean that a small change in the kite parameter destroys 
the orbit of (1/q, 1) when pq is odd, but not when pq is even. We prove the 
above statements in [S2]. 

My program Billiard King allows the reader to draw these pictures 
for any (smallish) rational parameter. Our monograph [S2] discusses these 
graphs in great detail. In fact, one could consider [S2] as an exploration of 
the arithmetic graphs associated to kites. 

2 You can download this program from my website http:/ /www. math. brown.edu/~res 
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3.2 The Regular Pentagon 



Let P be the regular pentagon. We scale P so that its vertices have the form 
i>j k for k = 0, 1,2,3,4. Here u = exp(2ni/5) is the usual primitive 5th root 
of unity. In this case, we have maps 7Ck : R J — > C given by 

4 
j=0 

Here X = (xo, x±, X2, X4). The map 7Ti is just the projection mentioned 
above. Note that iik and n^-k agree up to complex conjugation. Thus, the 
only remaining interesting projection is 7r 2 . 

Relative to any convex polygon, any periodic point x lies in a convex 
polygon P x consisting of points that have exactly the same dynamics. The 
outer billiards map moves the tile P x around isometrically. The pictures we 
show of outer billiards will draw the orbits of these tiles, rather than the 
orbits of individual points, because the pictures are more revealing. Figure 
3.3 shows the picture for the regular pentagon. 




Figure 3.3: Some of the periodic orbits for the regular pentagon 
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In the case of the regular pentagon, there are two kinds of orbits, those 
lying in regular pentagon tiles and those lying in regular 10-gon tiles. This 
is worked out in [T2], and also discussed in [BC]. The smaller the tile, the 
longer the periodic orbit. 

The left half of Figure 3.4 shows a plot of tc 2 (T(P,p)), where p is a point 
of S that lies inside a pentagonal tile. We choose the period to be fairly 
long, so that the actual size of the graph is quite large. One can see the 
emerging fractal structure in the rescaled picture. The right half of Figure 
3.4 shows the corresponding picture for a periodic point that lies inside a 10- 
gon tile. The reader can probably see that these tiles somehow fit together 
"hand-in-glove" . 




Figure 3.4: Some of the periodic orbits for the regular pentagon 



Any choice of long periodic orbit will yield a picture that is essentially 
identical to one of the two above. In the case of the regular pentagon, the 
arithmetic graphs lie within a thin tubular neighborhood of a 2-plane in R 5 . 
Any other projection of the graph would look like an afline image of one of the 
two pictures shown above. This is why we say that, in the regular pentagon 
case, there are two fractal curves that somehow control the structure of the 
pinwheel dynamics. 

As we remarked in the introduction, we do not offer proofs of these results, 
because we will give complete proofs of the analogous results for the regular 
octagon. 
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3.3 The Regular Octagon 

Let P be the regular octagon. We scale P as in §2.4[ and we define the 
projection maps as in §3.21 In the case, there are two interesting projections, 
7T 2 and 7r 3 . (It turns out that 7r 4 always maps the graph into a line segment.) 

As we discuss in the introduction, in connection with Figure 1.3, we call 
a periodic orbit odd if it intersects Ej, the half-strip from Figure 2.4., in 3 k 
points, for k odd. It turnso out that "half" the periodic orbits that start in 
Sq have this property. The odd periodic orbits are the lightly colored ones in 
Figure 3.5 that happen to intersect Ej. This includes all the lightly colored 
octagons that lie outside the first large ring of 8 octagons. 




Figure 3.5: The tiling associated to the octagon 
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Figure 3.6 shows some ir^ projections of the arithmetic graphs associated 
to odd orbits. For 7r 3 , the projections corresponding to the even orbits look 
similar. Note that Figure 1.2, which was drawn using the subdivision method 
discussed in §4.2[ emerges as we rescale the pictures. 




Figure 3.7: The snowflake emerges. 



Now we are going to show the 7r2 projections of the same orbits. 
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Figure 3.7: The carpet emerges in the odd case. 



For the even orbits, the 7r 2 projection is an open curve that looks fairly 
similar. We refer the reader to the program OctoMap 2, where one can see 
much better pictures, in both the odd and even cases. 

3.4 The Regular Heptagon 

We can make the same constructions for the regular heptagon as we made 
for the regular pentagon and octagon. The analogue of Figure 3.3 exists, 
but is not yet understood. The two interesting projections in this case are 
7r 2 and 7r 3 . Unlike the two cases we considered above, there is an explosion 
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of distinct pictures in this case. Sometimes the 7r 2 projection looks nicer 
and sometimes the 71% projection looks nicer. Sometimes both projections 
look incomprehensibly complex. Here are two simple examples, both 712 
projections. 




Figure 3.10: Some examples associated to the regular 7-gon. 



These examples do not even begin to capture the vast array of pictures 
one sees for the regular 7-gon. We encourage the experimentally-minded 
reader to explore the situation for themselves. 
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4 The Snowflake and the Carpet 



4.1 The Hausdorff Topology 

Let Ki,K2 C R 2 be two compact subsets. We define S(Ki,K2) to be the 
smallest e such that K\ is contained in the e-tubular neighborhood of K2 
and vice versa. The function 5 is known as the Hausdorff distance between 
Ki and K 2 . Given any compact subset VI C R 2 we let M(Q) denote the set 
of compact subsets of VI. Equipped with the function 8, the space M(Q) is 
itself a compact metric space. 

When we have a sequence of uniformly bounded compact subsets {K n } 
and we say that it converges, we mean that it converges in M(Q) for some 
large compact that contains all the individual sets. The set f2 just serves 
as a kind of container, so that we can speak about the convergence as taking 
place within a compact metric space. The notion of convergence we get 
is independent of the choice of fl Indeed, we could say more simply that 
K n ->■ K iff d(K n , K) -> 0. We call K the Hausdorff limit of {K n }. 

One case of interest to is is when we have a sequence {P n } of polygons 
whose diameter tends to 00. If we can a compact subset K and a sequence 
of similiarities {S n } such that S n (P n ) — > K in the Hausdorff topology, then 
we call K a rescaled limit of {P n }- In this case, the contraction factor of S n 
necessarily tends to 0. 

4.2 The Snowflake 

We call the snowflake J3. Figure 4.1 shows the subdivision of a right-angled 
isosceles triangle into 5 smaller ones. If the long side of the triangle has 
length 1, then the square on the right hand side of the figure has side length 
\f2 — l. The three small triangles in the middle have the same size. 




Figure 4.1: Subdivision rule 



To produce the snowflake, we start with the union of 8 isosceles triangles 
shown at the top right in Figure 4.2. Then we apply the subdivision rule 
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iteratively. The Hausdorff limit is our snowflake. We denote this limit by I 3 . 




Figure 4.2: Producing the snowflake by iteration. 



One nice feature of our construction is that 7 3 contains all the vertices 
of the triangles at every stage of the construction, and the union of these 
vertices is dense in I 3 . Moreover, all these vertices lie in Z[ou], where u is the 
usual primitive 8th root of unity. The triangles in our pictures are 2 colored 
in a natural way, according as to whether they point inward or outward. 
Assuming that we label a triange if it points inward and 1 if it points 
outward, our subdivision rule changes the colors according to the following 
rule: 

0^ 10001; 1 ^ OHIO. (10) 

We will see the significance of this structure below when we examine the 
carpet. 
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4.3 A Hidden Symmetry of the Snowflake 

For any right isosceles triangle T, let L(T) denote the curve one obtains by 
iterating the basic snowflake substitution rule and taking a limit. Assuming 
that T has a distinguished side a, let L(T, a) denote the portion of L(T) 
whose endpoints are the endpoints of a. 




Figure 4.3: Some special triangles 

Let T, Ti, and T 2 be the three right isosceles triangles shown in Figure 
4.3. The two triangles Ti and T 2 have the same size as each other, and are 
\pl — 1 times as large as T. Let A = L(T,a), where a is the edge of T 
bounded by the two white vertices. Let A\ = T(Ti,<7i), where o\ is edge of 
Ti bounded by the white and black vertex. Define A 2 = L(T 2 , a 2 ) similarly. 

Lemma 4.1 A = A 1 U A 2 . 

Proof: Let d be the Hausdorff distance between A qnd AiVJ A 2 . The 
substitution of T consists of 5 smaller triangles, T[, ...,T 5 ', with T[ = Ti and 
T' 2 U Tg related to T 2 just as Ti U T 2 is related to T. There are edges a' 2 and 
<7 3 of T 2 ' and Tg such that o 2 = a' 2 U a' 3 . The remaining triangles T' A and T^ 
lie to the left of the bottom white vertex of T. 

Since Ti = T[, we have A-y C A. At the same time, the pair of arcs 
(A,AiUA 2 ) is similar to the pair of arcs (A 2 , A' 2 U A' 3 ), where A' 2 = L(T 2) a' 2 ) 
and A' 3 = L(T 3 , a' 3 ). From this we see that the Hausdorff distance between A 
and A\ U A 2 is the same as the Hausdorff distance between A 2 and A 2 U A' 3 . 
But, by scaling, the latter distance is (\/2 — l)d. So, we have the equation 

d= (V2-l)d, 

which of course forces d — 0. 4 
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4.4 The Carpet 



We call the carpet Ii- The carpet is produced by substitution rule that is 
combinatorially identical to the one that produces the snowflake. This time, 
we have two shapes, parallelograms and trapezoids. The initial shapes have 
their vertices in Z 1 . Our figures below are accurate. 

Figure 4.4 illustrates how the parallelogram P is replaced by the sequence 

S P (P) =T 1 U? 2 U?3UP 4 UT 5 

of 5 smaller parallelograms and trapezoids. 




Figure 4.4: Substitution rule for the parallelogram 

Figure 4.5 shows how the trapezoid T is replaced by the sequence 

S T (T) = ? 1 UT 2 UT 3 UT 4 U? 5 

of smaller parallelograms and trapezoids. We say sequence here rather than 
union because the pieces in the sequence St{T) overlap each other. To make 
the substitution clear, we first draw Pi U T 2 U T 3 . Then we add P 4 and P 5 . 
The partition is invariant under reflection in the vertical line of symmetry. 




Figure 4.5: Substitution rule for the trapezoid 
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Let X be a cyclically ordered finite list of parallelograms and trapezoids, 
similar to the ones above. We let \X\ be the union of shapes in X. We let 
Sp(X) denote the result of substituting Sp(P) for all parallelograms P G X. 
Likewise we define St(X). Note that Sp(X) and St(X) both have natural 
cyclic orders, inherited from the ordering on X. Our carpet is 

I 2 = lim \I 2 {n)\ I 2 (n) = S P (S T (I 2 (0)). (11) 

Here ^(0) is a union of 8 copies of T arranged in a square pattern that winds 
around twice. (We use 8 rather than 4 so that the seeds for I 2 and J 3 have 
the same cardinality.) Figure 4.6 shows the sets I 2 {n) produced by iterated 
substitution, for n — 0, 1, 2, 3. Since the pieces in overlap, there is more than 
one partition we could draw. We have chosen the partition in which every 
trapezoid is drawn above every parallelogram. 




Figure 4.6: Iterating the substitution rule 
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4.5 Another View of the Carpet 

Figure 4.7 shows another substitution rule that generates J 2 - In Figure 4.7, 
we are showing how to replace a square by a subset of a square. Here we 
think of the dark triangles and quadrilaterals as markings that tell us how 
to orient the pieces when we iterate the basic rule. We obtained this rule 
by looking at Figure 4.6. The fractal obtained from this alternate rule is a 
subset of the Sierpinski carpet. It is obtained from the Sierpinski carpet by 
systematically deleting certain squares. 




Figure 4.7: An alternate substitution rule 



Lemma 4.2 I 2 is the fractal produced by the alternate substitution rule. 

Proof: We will start with the original rule and keep modifying it until we 
get to the alternate rule. Along the way, we show that the modifications 
don't change the limiting fractal. 
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For starters, we divide each parallelogram in half. In this way, we interpret 
the original substitution rule as a being generated by a rule for trapezoids 
and a rule for isosceles triangles. The left hand side of Figure 4.8 shows the 
new Sp(P), The right hand side of Figure 4.8 shows the new St(T). 




Figure 4.8: The first modification 



We have redrawn the second halves of Figures 4.4 and 4.5. In the case of 
Figure 4.5, we have drawn all the trapezoids on top. In drawing this piece, 
we arbitrarily chose to draw the left (divided) parallelogram on top of the 
right one. The substitution rule for the isosceles triangle does not respect 
the symmetry of the triangle. We record the asymmetry by orienting the 
hypotenuse of each triangle, as shown. The new substitution rule clearly 
generates J 2 . We still call our new rules Sr and St- Here R stands for a 
right isosceles triangle and T still stands for the trapezoid. 

Figure 4.9 shows two successive modifications of the rule for St- 




Figure 4.9: Two modifications of St- 



For the first modification, we add another trapezoid to St(T). This 
new trapezoid appears in Sr(St(T)), so the new rule rule generates the 
same limit. We give the same names to our new substitution rules. For the 
second modification, we delete two of the isosceles triangles. Call the deleted 
triangles R\ and R 2 . Observe that Sp(Rk) C S T (T) for k = 1,2. Therefore, 
the latest subdivision rule has the same limit as the original. We again call 
the generators of the rule St and Sr. This time St and Sr replace one piece 
by a union of pieces that have pairwise disjoint interiors. 

From here it is a simple exercise to check that Sr o St implements the 
rules shown in Figure 4.7. 6 
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4.6 The Map from the Snowflake to the Carpet 

The initial seed for I 3 is a cyclically ordered list h(0) of 8 isosceles triangles 
Ti, r 8 . The initial seed for I 2 is a cyclically ordered list h(0) of 8 trapezoids 
Ti, T 8 , which wrap around twice, as we have already mentioned. We have 
a correspondence Tj -h- Tj between the two seeds. 

Since the two substitution rules are combinatorially identical, we induc- 
tively get a surjective map from the triangles in h{n) to the quadrilaterals in 
I 2 (n). This map respects the cyclic ordering on both sets, and it also respects 
the relation of descendence. Say that a shape r n of Ij(n) is a descendent of 
a shape r m of Ij{m) if the substitution rule, applied n — m times, to r m , 
produces a list of shapes that contains r n . 

We define : J 3 — > I 2 by taking the limit of our correspondence. Given 
x G I3 we can find a sequence {r n } of triangles, with r n G X 3 (n), that accu- 
mulates on x. We define 0(x) to be the limit of the corresponding quadrilat- 
erals in I 2 . 

Lemma 4.3 is well-defined. 

Proof: Suppose that {r' n } is another sequence of triangles that converges to 
x G 1 2 . Let S n and S' n be the shapes of I 2 (n) corresponding respectively to 
r n and r' n . It suffices to prove that diam(S'„ U S' n ) — > 0. 

Call two triangles in Izirn) close if they are either identical or else con- 
secutive in the cyclic order. Let d n be the largest integer such that r n and 
T' n are descendents of close triangles in h(d n ). Since the distance between r„ 
and r' n converges to 0, and 1% is evidently an embedded curve, d n — > oo. 

We make the same definitions as above for I 2 . We can find close shapes S n 
and T*' n in I 2 (d n ) such that S n is a descendent of S n and S' n is a descendent of 
Ti' n . Note that S n and have intersecting boundaries in all cases. Also, the 
diameters of these shapes tends to with n. Hence diam(S Tt U S'J — > 0. Fi- 
nally, we have S n C S n and S' n C T! n . Hence diam(S' n US'' t ) — > as desired. 4ft 

Essentially the same proof shows that is continuous. Since all the ap- 
proximating correspondences are surjective, is surjective as well. 

Remark: Note that we could get an equally canonical map if we compose 
with some isometry of I 2 or 1%. This basic (though trivial) ambiguity comes 
up in Statement 3 of our Main Theorem below. 
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4.7 Geometry of the Map 

Now we give some geometrical information about the map (f) : ^3 — > I 2 defined 
in the previous section. Note that I3 contains the vertices of triangles in I^{n) 
for all n. Say that a special point of ^3 is a right-angled vertex of I^{n) for 
some n. The set A3 of special points is dense in I 3 . 

Say that a special point of I 2 is a midpoint of the long edge of some 
trapezoid in I 2 (n). We could equally define the special points to be the 
centers of the parallelograms. The two definitions coincide. Let A2 denote 
the set of special points of I<i- The set A2 is the set of midpoints of the 
edges of the countable family of squares that arises from the substitution 
rule explained in §4.51 

Lemma 4.4 <f>(A 3 ) = A 2 . 

Proof: Let v G A3 be some point. Then v is the right-angled vertex of 
a triangle tq of Isfo) for some n. Note that v is also the right vertex of a 
triangle r k of I^{n + k) for k = 1,2,3... The triangles ti,t 2 , ... all have the 
same type, either inward pointing or outward pointing. 

Consider the outward pointing case first. In this case the quadrilateral 
T k of 12(11 + k) corresponding to r k is a trapezoid. Inspecting our subdivision 
rule, we see that T k+ i is the middle trapesoid of 5r(7fc). Looking at Figure 
4.5, we see that T k and T k+ i have a common distinguished point. This holds 
for all j. Hence f]T k = <p(v) is the common distinguished point of all these 
trapezoids. Hence 4>{v) G A 2 . 

In the inward pointing case, the quadrilateral P k corresponding to r k is a 
parallelogram. Here P k +i is the middle parallelogram of Sp(P k ). But these 
parallelograms have a common center. The intersection f]P k = <p(v) is this 
common center. Looking at Figure 4.5 again, we see that the centers of the 
parallelograms coincide with distinguished points of trapezoids. So, again 
(j)(v) G A 2 . 

Our argument so far shows that <j)(A^) C A%. Any x G A2 is the distin- 
guished point of some trapezoid T of some ^(n). But then there is a triangle 
r of Iz(n) which corresponds to r. By construction x = <f)(v), for the right 
vertex v of r. Therefore 4>{A 3 ) = A 2 . 4* 

Let B 3 C I3 denote the set of points that are acute vertices of triangles 
in I$(n) for all the different n. Then B3 is dense in I3. Let B>2 C I2 denote 
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the set of vertices of quadrilaterals in I 2 (n) all the different n. The set B 2 is 
the set of corners of the distinguished squares mentioned above in connection 
with A 2 - The curve I3 looks "oscillatory" in the neighborhood of any 1 G B3. 
Any line through x intersects ^3 infinitely often in every neighborhood of x. 

Lemma 4.5 (f>(B 3 ) = B 2 . 

Proof: The proof is similar to the one for the previous result, so we will be 
a bit sketchy. Any x G B 3 (n) is the acute vertex of a triangle in J 3 (n + k). 
The triangle t^+i is either the first or last triangle in the list of 5 triangles in 
I 3 (n + k + l) that replaces rfc. Whether r^+i is "first" or "last" is independent 
of k. 

Let Tfc be the quadrilateral of I 2 {n + k) corresponding to rfc. By construc- 
tion, these quadrilaterals all share a common vertex, and <fr(x) must be this 
common vertex. This shows that <p(B 3 ) C B 2 . The proof that B 2 C <p(B 3 ) is 
just as in the previous result. 4 

Remark: It is worth pointing out one more feature of the map 0: It is far 
from 1 to 1. The set I 2 {n), considered as a list of quadrilaterals, contains 
multiple copies of the same trapezoid. Many trapezoids appear 2 n times in 
I 2 (n). From this we see that, for any n, there are 2 n segments of ^3 that cj) 
identifies. For this reason, the canonical surjection from ^3 to I 2 is a many-to- 
one map. We think if it is a kind of fractal version of the universal covering 
map from the line to the circle. 

4.8 A Hidden Symmetry of the Carpet 

The snowflake has a hidden symmetry, as discussed in §4.31 Given the close 
correspondence between the carpet and the snowflake, we would expect a 
combinatorially identical hidden symmetry to appear for the carpet. This is 
indeed the case. In this section we describe the symmetry and sketch the 
proof. 

The hidden symmetry of the carpet manifests itself in two ways, and this 
at first seems different from what happens with the snowflake. However, in 
the case of the snowflake, there are really two kinds of right-angled isosceles 
triangles, so the symmetry there actually does arise in two ways. 
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Figure 4.7: Some special quadrilaterals. 



The left side of Figure 4.7 shows a large parallelogram and two smaller 
shaded trapezoids. We can generate a subset of a suitably scaled copy of the 
carpet starting with the large parallelogram as a seed and then considering 
the portion of the resulting limit that connects the two white dots in Figure 
4.7. Let A denote the resulting set. Likewise, we start with each shaded 
trapezoid as a seed, and consider the portion of the resulting limit that 
connects a black dot to a white dot. Let A 1 and A 2 be the resulting sets. 

The right hand side of Figure 4.7 shows a trapezoid and two shaded 
parallelograms. We let B and B\ and B 2 be the portions of carpets produced 
by the same construction as we did for the As, but interchanging the roles 
of black and white. 

Lemma 4.6 A = A x U A 2 and B — B 1 U B 2 . 

Proof: The proof here is formally isomorphic to what we did for the snowfiake 
We omit the details. 6 
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5 The Main Result 



5.1 Statement of the Result 

Let lo = exp(27ri/8) be the usual primitive 8th root of unity. Let P be the 
regular octagon whose vertices are powers of u. Let Eo be the pinwheel strip 
associated to P as in §3.31 Let Eq C E denote the half-strip that lies to the 
right of the line x = 1 + y/2. See Figure 5.1. 




Figure 5.1: The relevant sets 



§7.51 we prove 

Lemma 5.1 Relative to the the pinwheel map for the regular octagon, any 
periodic orbit that intersects Eq nontrivially does so in exactly 3 k points for 
some k = 0, 1, 2.... 

Given a periodic point x e Eq, we define \x\ to be the exponent k such 
that the pinwheel orbit of x intersects Eq exactly 3 fc times. We call a sequence 
{x n } G Eq a good sequence if \x n \ — > oo as n — » oo and he congruence of \x n \ 
mod 4 is odd, and independent of n. 



37 



Recall that 7r 2 and 7r 3 are the projections from R 8 to R 2 defined in the 
previous chapter. Define 

s 2 = V3; s 3 = l + V2. (12) 

Let : J3 — > I 2 be the canonical surjection described above. 

Theorem 5.2 (Main Theorem) Let {x n } be a good sequence. Also, let 
A n C R 8 be the arithmetic graph of x n . Then the following is true. 

1. 7i 2 (A n ) and 7r 3 (A n ) are closed polygons for all n. 

2. Let Sk-n be dilation by about the origin. Then the following limit 
exists. 

T k = lim S kn o ii k (A n ). 

n— >oo 

3. r 3 = I 3 , where J3 is scaled so that one of the 8 isosceles triangles in its 
seed has vertices 

(0,0); (l/2,-s 3 /2) (-s 3 /2,-l/2). 

4. T 2 = h, where I 2 is scaled so that its center is (3/2,3/2) and one of 
its main corners is (—3/2,3/2). 

5. Identifying points on ir s (A n ) and TT 2 (A n ) if they come from the same 
point on A n , we get a polygonal surjectin 

4> n ■ S 3tn ir 3 (A n ) S 2in Tl 2 (A n ). 

We have lim^oo 4> n = F o (f>. Here F is an isometry of I 2 and <fi : I3 — > 
I 2 is the canonical surjection from the snowflake to the carpet. 

Remarks: 

(i) Experimentally, we see that Theorem 15.21 also holds when we consider se- 
quences of points in S — Sj, but these cases lead to annoying complications. 
We are interested in establishing a nice robust version of the phenomenon, 
but not necessarily the sharpest possible version. 

(ii) The limits are exactly the same when \x n \ = 1 mod 4 and when \x n \ = 3 
mod 4. The only thing that changes is the global isometry F. In the former 
case, F is orientation preserving (and can be taken to be the identity) and 
in the other case F is orientation reversing. 

The rest of the paper is devoted to proving Theorem 15.21 
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5.2 Restricting the Domain 

Let P be the regular octagon, and let S be the pinwheel strip described in 
§2.41 Consider the octagons 

Of = P + A;(v / 2 + 1,±1); £ = 1,2,3... (13) 

Of and 0^7 share a common vertex and fit inside So as shown in Figure 
5.2. Figure 5.2 shows the cases k = 1,2, but the remaining cases look the 
same. 




Figure 5.2: The region R\. 



We let Rk be the region of S bounded by Of U 0^7 and Of +1 U 0^ +1 . 
Figure 5.2 shows the region Ri, which is the main region of interest to us. 
Note that Rk is translation equivalent to R\ for all k > 1, so the picture in 
Figure 5.2 is typical. 

The first purpose of this chapter is to prove the Invariance Lemma. This 
result implies that each Rk is forward invariant under the pinwheel map. 
Following this, we will prove 

Lemma 5.3 Each Rk is an invariant set for the pinwheel map. Furthermore, 
if Theorem \5.S\ is true for sequences in R\ then it is true as originally stated. 
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5.3 The Pattern of Octagons 

let P be as above. The octagons we described above are part of the infinite 
pattern of octagons shown in Figure 5.3. Figure 5.3 shows a part of an 
infinite set of octagons that, as it turns out, is invariant under the outer 
billiards map. Figure 5.3 also shows the strip E . 




Figure 5.3: Octagonal orbits 

The 8 dark rows of octagons are obtained by translating the central one 
by a vector of the form 

ku n (V2 +1,1); k,ne Z; to = exp(2ni/S). (14) 

The remaining octagons interpolate between these dark ones. The octagons 
are situated in such a way that the outer billiards map is entirely defined on 
the interior of each one. The octagons come in an infinite family of rings, or 
necklaces N 1: N 2: N 3 , .... The octagons are precisely the two octagons of 
S n N k . Compare [GS] and [K]. 
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Lemma 5.4 Each necklace is an outer billiards orbit. Moreover, the region 
between consecutive necklaces is invariant under the outer billiards map. 

Proof: We change notation slightly. Let 0(n, k; 0) be the octagon we get by 
translating the central one by the vector in Equation[T31 Let 0(n, k; i) denote 
the octagon that is i "clicks" away from O(n,k;0) going counterclockwise 
around iV fc . 

The outer billiard map has the same action on all the octagons lying on 
the vertical line segment between O(—l;k;0) and 0(l,k; — 1) because such 
octagons all lie in the sector shown in Figure 5.4. The outer billiards map 
reflects each of the octagons under construction through the apex of the 
sector in Figure 5.4. We check easily that the action is as follows. 

O(-l,fc;0) -> 0(4,fc;l); 0(1, k; -1) ->■ 0(5, k; 0). (15) 




Figure 5.4: Octagons in a sector 

The octagons 0(4, k;l) and 0(5, k;0) have their centers on the same 
vertical line segment L, as shown for k = 2 in Figure 5.4. By symmetry, the 
octagons between 0(— 1, fc; 0) and 0(1, k; — 1) map to the octagons of the 
grid that have their centers on L. These octagons all belong to Nk. For the 
remaining octagons in Nk, the result follows from the 8-fold symmetry of the 
picture. 

Call two octagons Oi G N k and 2 G N k+ i close if they are as close 
together as possible. Two close octagons determine a polygonal region "be- 
tween them" like the one that is shaded in Figure 5.4. Call such a region a 
buffer. If the two octagons belong to the sector from Figure 5.4, then the 
buffer is mapped between N k and N k+1 . But then, by symmetry, all buffers 
are mapped between N k and N k+1 . But every point between N k and N k+ i is 
contained in a buffer. 4k 
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5.4 Invariant Domains in the Strip 

Let E denote those points in So that lie to the central octagon P. A point 
p lies in E if every vector pointing from P to p has positive x coordinate. 
Let $' : E — >■ E be the first return map of the <ft 2 , the square outer billiards 
map. 

Corollary 5.5 For k > 2 the set R k is invariant under 

Proof: For k > 2, the region between A^ and Afc +1 intersects E exactly in 
the set Rk- 4 

Remark: The region between Ni and N2 intersects Eo in a set that is some- 
what larger than Ri., as one can see by looking carefully at Figure 5.3. It is 
for this reason that Corollary 15.51 fails for Ri. 

Corollary 5.6 For k > 4 the set Rk is invariant under <E> ; the pinwheel map. 
Proof: For k > 4 the argument in Lemma [27T1 goes through, and shows that 

Below we will prove the following result. 
Lemma 5.7 Let fx = 1, 2, 3... Lei p E Ri and 

q = p + 2fu(V2+ 1,1) G i? 1+A . 
T/ien i/ie arithmetic graphs T{p) and coincide. 
Corollary 5.8 The set R^ is ^-invariant for all k > 1. 

Proof: By Lemma [5771 translation by the vector 4(v2 + 1, 1) conjugates the 
action of $ on Ri to the action of $ on R 5 . We already know that R$ is 
invariant under $ and the conjugation property implies the same result for 
Ri- A similar argument takes care of R2 and R3. We have already handled 
the remaining cases. 4k 
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Remark: Note that $ and $' do not coincide on R\, because Corollary 15.51 
is false for Ri whereas Corollary 15.81 is true. These maps in fact coincide for 
k = 2,3,4.... 

Proof of the Invariance Lemma: We want to show that E($ is $ invariant. 
Let N£ be the set of half-octagons obtained by taking the right halves of all 
the octagons in Since <ft 2 is a piecewise translation, the set N£ is invariant 
under <fi 2 . But then the set 

s k = N+ n s 

is $' invariant for k > 4. Hence, by the argument in Lemma 12. 1[ the set Sk 
is ^-invariant for k > 4. As in the proof of Corollary 15.81 we now conclude 
that Sk is $ invariant for all k = 1, 2, 3... Finally 

Sq = Sk U Rt+i U Sk+i U Rk+2 
We have decomposed Eq into $ invariant sets. Hence Eq is also <3> invariant. 

Proof of Lemma l5.3t Suppose we want to prove Theorem l5.2l for a sequence 
of periodic points {x n } in Eq. Since $ is the identity on each of the orbits 
Nk, we have x n G Rk n for large n. Chopping off the initial portion of the 
sequence, we can assume that x n G Rk„ for all n. 

By Lemma 15.71 we can find a point y n G i?i U i?2 such that x n and y n 
have the same arithmetic graph. Hence, it suffices to prove Theorem 15.21 for 
the sequence {y n } G R\ U i?2- 

The intersection i?i U R2 is a sing le point, namely C = (2 + 2y/2, 1). We 
call two points specially related if reflection through ( interchanges the two 
points. We check by hand the following property, which we call Property P: 
If y and z are specially related then $(?/) and $(2) are specially related and 
+ 6(z) = 0. 

From Property P we see that the arithmetic graphs of y and z are iso- 
metric to each other. Using this basic principle, we can replace each y n G R2 
by the specially related point z n without changing the isometry type of the 
rescaled limits of the arithmetic graphs. 4k 

Remark: In §[7] we will explain in detail how we actually check Property P. 
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5.5 Equivalent Points in Strips 



The rest of the chapter is devoted to the proof of Lemma 15.71 In this section 
we make a general observation about strip maps. In the section following 
this one, we will apply our observation to the case of the regular octagon. 

Let S x and S 2 be non-parallel strips. Let V and Wi,W 2 be the three 
vectors shown in Figure 5.3. We say that p, q G Si are related if p — q = kW\ 
for some integer k. Likewise, we say that p, q G S 2 are related if p — q = kW 2 
for some integer k. 

c 



B 



A 




Figure 5.4: A Strip map 

Lemma 5.9 Let f be the strip map associated to (S 2 ,K). Ifp,q G Si are 
related then f{p),f{q) G S 2 are related. 



Proof: Our statement is an affinely invariant one. So, we normalize our 
strips so that they are perpendicular and each have width 1. In this case, / 
translates the horizontal squares A,B,C, ... to the vertical squares with the 
same labels. If p, q G Si are related then they are in the same positions rel- 
ative to the horizontal squares that contain them. But then f(p), f(q) G S 2 
are in the same position relative to the vertical squares that contain them. 
Hence f{p) and f(q) are related. 4 
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5.6 Proof of Lemma 15.7 

We first prove Lemma [5.71 in case that A = 2/i is even. 
We define 

W 2n = W 2n +i = 2u n (V2 + 1,1). (16) 

The vector W = W\ is parallel to both pointed strips S and Si. The vector 
W 2 = W 3 is parallel to both S 2 and S 3 , and so on. 
We say that two points p, q G T>k-i are equivalent 

q-p = fiW k -i\ (JteZ. (17) 

We let fik-i(p,q) be the integer in Equation [T71 

Let fk be the strip map associated to V k ), as in §2.31 Let rrikip) be 
as in Equation [2j That is, 

fk{p) = p + m k {p)V k . (18) 

Lemma 5.10 Let f = f k . Ifp,q G S fc _i are equivalent and k is odd, then 
f{q),f{q) are equivalent in S& and 

M/(p), /(?)) = ^*-i(p> ?); m fcb) = ™*(?)- 

Proof: By rotational symmetry, it suffices to consider the case when k = 1. 
Suppose that k — 1. Note that S fl Si is a strip parallel to Wo = W\. So, 
p G S fl Si if and only if q G S fl Si. Suppose that p,q E So fl Si. Then 
fip) — V an d f(q) = q- Also, Wi = Wo- Hence f(p) and f(q) are equivalent 
in Si and 

MpI /(?)) = //o(p, /i(p) = = o. 

If p, q G S — Si then p + V\ G Si and likewise q + V\ G S x . Here 

A(?) - AGO = q-p- 

Hence 

A*i(/(p)> /(?)) = A*o(p, 9); /i(p) = A(?) = L 

4 
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Lemma 5.11 Let f = f k . Ifp,q € £&_! are equivalent and k is even, then 
f(q),f(q) are equivalent in and 

Vk(f(p), f(q)) = Pk-i(p, q); m k(q) = m k {p) + n k -i(p, q)- 

Proof: By rotational symmetry, it suffices to consider the case k = 2. In 
this case, the strips Si and £2, as well as the vectors W\ and W2 and V = V2 
are related precisely as discussed in §5.51 Lemma 15.91 now implies that p and 
q are equivalent in S 2 an d ^{Pi q) — Hi{p-,q)- 

Referring to Figure 5.4, which is an affine image of the octagon picture, 
the number fix (p, q) counts the number of parallelograms one needs to shift 
in order to move p over to q. For instance, if p were in parallelogram A and 
q were in parallelogram C in Figure 5.4, then Ho{p, q) = 2. From Figure 5.4, 
we see clearly that m 2 (g) = m k {p) + p,i(p, q). 4 

Now we can finish the proof of Lemma 15.71 Suppose that p, q G So are 
equivalent. Let \i = fio{p,q). Then, referring to Equation HI we have 

8 

$(g) - $(p) = £(m fc (g) - m h (p))V k = l ^ fiV k = 2 0. (19) 

i=l k even 

Equality 1 comes from Lemmas 15.101 and I5.11[ and Equality 2 comes from 
the fact that Vq— = V 2 and V$ = — V4. 
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6 A Toy Model 



6.1 The Map 

The work in the previous chapter reduces the proof of Theorem 15.21 to the 
study of the pinwheel map $ : R\ — > Ri, where R± is shown in Figure 5.2. 
This first return map turns out to be a polygon exchange map. In this 
chapter, we will briefly discuss a simpler but related polygon exchange map. 

Though the simpler map is not the one we will ultimately study in our 
proof of Theorem I5.2[ we think that it will prepare the reader for the more 
complicated system we do study. First, the basic shapes that arise in the 
simple system here arise in the more complicated system we study later. 
Second, both the system here and the one we do study have a renormalization 
scheme. It is the renormalization scheme for the first return map that leads 
to the self-similar nature of the arithmetic graphs. 

The dynamical system we study here is well known, and indeed also arises 
in the study of outer billiards: It is the first return map to a certain invariant 
domain for a dynamical system generated by the outer billiards map and the 
rotational symmetry group of the regular octagon. Compare [T2] and [BC]. 
The system is a self map of the kite-shaped region X, shown as the shaded 
region in Figure 6.1. The region X is best defined in terms of a large regular 
octagon. The smaller regular octagon, which is a subset of X, is drawn in 
for later reference. 



Figure 6.1: The region X. 

As shown in Figure 6.2, the region X can be partitioned in two ways into 
right angled isosceles triangles, 

X = Ax U B 1 = A 2 U B 2 . (20) 
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There is a unique orientation-preserving isometry '■ Ai — > A 2 . The 
map 4>a rotates 37r/4 radians clockwise about the center of the small octagon 
in Figure 6.1. Likewise, there is a unique orientation preserving isometry 
4>b '■ B\ — > B 2 . The map 4>b rotates 7r/4 radians counterclockwise about the 
center of the big octagon in Figure 6.1. We have a map : X — > X which 
restricts to (ft a on A 1 and <p B on B ± . The map is defined except on the 
segment common to A 1 and B x . 




Figure 6.2: Two partitions of X. 



6.2 The Renormalization Scheme 

Let O C X denote the inner octagon in Figure 6.1. By construction, the 
center of O is a fixed point of (p, and every point of O has period 8 with 
respect to 0. 




Figure 6.3: Renormalization Scheme 
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Let X' be the shaded component of X — O shown in Figure 6.3. We have 
partitioned X' into two pieces, which we call A[ and B[. The third power 3 
maps X' to itself. Figure 6.2 shows the pieces A 2 = 0(^4'i) and A' 3 = <p{A 2 ) 
and B' 2 = <P(B[) and B' 3 = <j)(B' 2 ). 

We have not drawn A' A ,B' 4 C X', but we observe that A' 4 and B' 4 relate 
to A[ and B[ in the same say that A 2 and B 2 relate to A\ and B\. More 
precisely, there is an orientation reversing dilation : X' — Y X such that 

0- 1 o0o0 = 3 (21) 

on X' . The map is the remormalization map. Having by itself is nice, 
but we have a second property that really pins things down. Namely, 

X-O = X'U0(X')U0 2 (X'). (22) 

In other words, every point of X — O lies in the orbit of a point in X' . 

6.3 The Periodic Points 

Our remormalization scheme allows us to get a complete understanding of 
the periodic points of <fi. The octagon O consists of points having period 8. 
(The center is the one point of O that has period 1.) For our purposes, it is 
nicer to think of O as having period 1 as a tile. The set O is preserved by 
<p. We set O = O. 

By equation I2TJ the region X' contains a small octagon 

O' = e-\O ) (23) 

that has tile-period 3. We set 

0' k+ i = <P(0' k ). (24) 

The tiles 0' k for k — 0, 1, 2 all have tile-period 3. They comprise a tile-orbit. 
Now we iterate. By Equation (2TJ the three octagons 

C^ = © _1 (0*); fc = 0,l,2 (25) 
all have tile-period 9. Setting 

O'U = #00, (26) 

we produce 9 small octagons Oq, O'l that all have tile-period 9. Next, we 

(3) 

produce 27 small octagons O k for k = 0, ...,26, all having tile-period 27. 
Figure 6.4 shows some of these octagons. 



49 



Figure 6.4: Cascade of periodic tiles 



Thus, we have produced, for each n — 0, 1, 2, 3..., a periodic octagonal tile 
having tile-period 3 n . We claim that every periodic point is contained in one 
of these octagons. To see this, we use Equation [22j Clearly, the only periodic 
points of order 1 lie in Oq. Suppose N is the smallest positive integer for 
which we have not proven our result. Let y' be a periodic point of period N. 

Since N > 1, we have y' G X — O. By Equation [22j the point y' lies in 
the orbit of some point x' G X'. Let x = 0(x'). Note that <f){X') and <f) 2 {X') 
are disjoint from X' . Since <f) N (x') = x', we must therefore have N = 3M for 
some integer M. But then, by Equation [2U we see that x has period M. By 
induction, x lies in one of our octagonal tiles. But then, by construction, so 
does x'. Now we know that x' is one of the periodic points we have already 
studied. So, the same goes for y', the original point of interest to us. 

Using the renormalization scheme, we have been able to classify all the 
periodic points of the system. Incidentally, we note that the complement of 
the octagons has dimension larger than 1, and the set of points where some 
iterate of is undefined has dimension 1. Therefore, has some aperiodic 
points. This is the same argument Tabachnikov [T2] gives for the regular 
pentagon. 
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7 The Pinwheel Dynamics 



7.1 The Partition 

We will work with the region Ri from Lemma [5.31 Both sides of Figure 7.1 
show R\. Note that R\ is the union of 4 copies of the region X from the 
previous chapter. 




Figure 7.1: The region inside the strip So- 



Figure 7.1 shows the action of $ on R\. The left hand side of the figure 
shows the partition of R\ into the maximal pieces on which $ and $ are 
constant, and the right hand side shows the images of these pieces under $. 
We will give coordinates for these tiles below. 

Just knowing that $ is a piecewise translation, we can nearly pin down 
the action of $ from Figure 7.1. The problem is that there are some pairs 



51 



of triangles, having identical shading, that are translates of each other. Such 
pairs are specially situated: One triangle in each such pair lies on the left 
hand side of R\ and the other triangle lies in the right hand side. To pin 
down the action of $ exactly, we mention that <3> maps any triangle on the 
left (respectively right) hand side of R\ to a triangle on the right (respec- 
tively left) hand side of R\. 

Remark: It seems worth also mentioning that $ maps any quadrilateral or 
pentagon on the left (respectively right) hand side of R\ to another quadri- 
lateral or pentagon on the left (respectively right) hand side of R\. Finally, 
$ fixes each of the central octagons. 

$ has one obvious symmetry and one hidden symmetry. Let p\ be the order-2 
rotation about the center point of R\. Then $ commutes with p. 

To describe the hidden symmetry, we observe that it is better to think of 
R\ as a subset of an infinite cylinder. Define a(x, y) = (x, y+2). The squared- 
map a 2 maps the bottom boundary of S to the top one. The quotient, 
obtained by identifying the two boundary components, is a cylinder. We 
think of Ri as a subset of this cylinder, so that certain pieces on the top of 
R\ are contiguous with pieces on the bottom. 

This point of view makes the picture look more symmetric. Consider 
the magenta pentagon K' shown on the middle left hand side of Figure 7.1. 
Let K" be the magenta triangle such that K — K' U K" is a kite, isometric 
to the rest of the kites in the picture. Q{K') lies at the very top of R\ 
and &(K") lies at the very bottom. Under the identification we have been 
discussing, &(K") is just the continuation of <&(K'). So, when we think of 
Ri as a subset of a cylinder, K is a single tile on which <3> is constant. The 
other pentagon/triangle pair on the left hand side of Figure 7.1 has the same 
analysis. Considered this way, we see that the partition of R\ consists just 
of kites, triangles, and octagons. 

Finally, we mention the hidden symmetry. Once we interpret R\ as being 
a subset of the cylinder, <£> commutes with a. 

Remark: It turns out that $ does not share the same symmetries. That is, 
it is not true that the maps $ and $ o p and $ o a coincide. However, as we 
will see, there is some relation between these maps. 
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7.2 Notation and Coordinates 



This section does not contribute to the mathematics at all. We include it for 
readers who would like to see explicit coordinates to that they can reproduce 
our experiments and calculations on their own. Let 

s = V2 + l. (27) 

Let O3 be the red octagon in Figure 7.1. The octagon 3 has radius 1/s. 
By this we mean that the distance from the center of 3 to any vertex is 
1/s. We point out another octagon, 4 , of radius 1/s 2 . The octagon 4 is 
homothetic to 3 , and the rightmost vertex of 4 coincides with the leftmost 
vertex of 3 . Both 3 and 4 are periodic tiles for $. See Figure 7.2 below. 
Now we distinguish 4 special points. 

• The point c\ = (3s, l)/2 is the center of R±. 

• The point C2 = (3s, 3)/2 is the center of the top half of R±. 

• The point c 3 = (2 + \/2, 2) is the center of 3 . 

• The point c 4 = (2\/2, 2) is the center of 4 . 

Let pk denote counterclockwise rotation by 2nk/8 about Ck- 




Figure 7.2: Some polygons of interest. 
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Figure 7.2 shows 4 shaded polygons. For these polygons, we will use a 
second labelling scheme in which O4 = B . Figure 7.2 shows the B labels. 
The union B2 U B 3 is the kite that lies to the left of the red octagon in Figure 
7.1. We have tried to shade B 3 in such a way that we point out O4 C B 3 . 

The vertices of 2B 3 are 

(5s-6,2s-l) (4s-4,4s-6) (3s-2,2s-l) (s + 3,3s-3) (2s + l,s + 2). 

(28) 

We divide these coordinates in half to get the coordinates of B 3 . 
The vertices of 2£?i are given by 

(2s, 4); (3 + s,3s-3); (3s-2,2s-l). (29) 

We divide these coordinates in half to get the coordinates of B 3 . 

We can deduce the remaining coordinates just from what we have already 
given. For instance, 

B A = p\{B 1 )- B 2 = pl(B 1 ). (30) 

The remaining polygons in the partition of R\ all have the form 

PT op?opf(B); ai ,a 2 e{0,l}; a 3 G {3,5}. (31) 

Here either B = Bj or B = B 2 U B 3 . 

Now we explain briefly how we verify that the return map $ is as stated. 
We explain this for B 3 . The check is the same for the remaining files. Recall 
that $ = /i6°"-°/i, where fk is the strip map associated to the pair 14). 
We verify that there are integers mi, ...,mi 6 such that 

• B\ := Bs + m^ G 

• B\ := B\ + m 2 V 2 G S 2 . 

• Bl := Bl + m 3 V 3 G S 3 . 

And so on. By convexity we just have to check this on the vertices of B 3 . 
This is a short calculation we omit. This calculation shows not only that 
$(£> 3 ) = _Bg 6 , but also that $ is constant on the interior of B 3 . 
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7.3 The Compressed System 

Let R C R\ denote the top half of R\. The region R is colored red and 
magenta on the left half of Figure 7.1. Note that R is a fundamental domain 
for the action of a on R\. We define a new dynamical system \l/ : R — > R as 
follows: ^(p) = if $(j>) e R and = a o if ^ i?. In other 
words, we use the action of o to move everything into R. 

Figure 7.3: The compressed system 

The compressed system has 22 regions. The regions Bq,...,B^ are as in 
Figure 7.2. For the record, here are the remaining regions. 

• B 5 = pKB,) and B 6 = p§(B 2 U B 3 ) and B 7 = p 5 3 (B 4 ). 

• B s = &{B X ) and B 9 = p 2 3 (B 2 U B 3 ) and B 10 = pl(B 4 ). 

• B k+ll = p\{B k )- A: = 0, ...,10. 

The regions B 5 ,B e , B-j all lie in the upper right corner of R. The regions 
B 8 ,B 9 , B w lie in the middle of R on the left side. The remaining 11 regions 
are symmetrically placed with respect to the first 11. Again for the record, 
the red colored regions are 

B i+llj ; i = 2, 4, 6, 7, 8, 9; J = 0,1. (32) 

We will use this notation system in the next chapter, just so that we can list 
things in a completely explicit way. 
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Figure 7.3 does for \l> what Figure 7.1 does for $. The left hand side 
shows the partition of R into maximal regions on which \1/ is constant. The 
blue polygons are the ones on which ^ = $ and the red polygons are the 
ones on which ^ = a o $. The right hand side shows the images of these 
pieces under \1>. Once we mention that \1/ is a piecewise translation that maps 
all the triangles on the left (respectively right) side of R to triangles on the 
right (respectively left) side of R, the partitions alone determine We also 
remark that ^ maps the kites and pentagons on the left (respectively right) 
side of R to the left (respectively right) side of R. 

The map ^ has a symmetry that it inherits from $, and also a new one. 
First, \l/ commutes with the rotation p 2 about the center point of R. Second, 
the partition on the right hand side of Figure 7.3 is a mirror image of the 
one on the left. (The line of symmetry joins the centers of the red octagons.) 
Finally, \1/ moves the one partition to its mirror image by way of translations. 

7.4 Classification of Periodic Orbits 

Lemma 7.1 p G R is a periodic point for ^ if and only if p is a periodic 
point for $ . 

Proof: We write 00 = ^0= P- Let <p n be the nth iterate of 0o under <£> and 
let ip n be the nth iterate of ipo under \1>. Since $ commutes with a, we have 

&» = <7**(^ n ) (33) 

where d n counts the number of k e {0, (n — 1)} for which ^ lies in a blue 
region in the partition of R (on the left side of Figure 7.3.) 

Suppose that p has $-period n. Then <p n = <f) . This means that 
o~ dn {jp n ) = V'o- Since both ip n and ipo belong to R, we must have d n even. 
Hence ip n = This shows that the p is periodic for \1>. On the other hand, 
suppose that p has \l/-period n. Then %p n = ip$. If d n is even then <p n = <J)q. 
If d n is odd, then <i 2 n = 2d n because ip n +k — ipk for each k — 0, ...,n. Hence 
02n = 0o- Hence p is ^-periodic. 4 

Figure 7.4 shows the beginning stages of a packing of octagons in R. As 
we will see momentarily, every periodic point of \& lies in one of the octagons 
in this packing. 
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Figure 7.4: Image of the central portion 



We will establish the result about the packing using a renormalization 
scheme that is very similar to the one used in the previous chapter. The 
remormalization scheme is based on the dilation 

Q(z) = az + b; a = (s + 1) + i(s + 1); b = -3i(s + 1) (34) 

We have Q(S) = R, where S C R is the portion of R between the two 
octagons 0\,Oi C R. The left hand side of Figure 7.4 below shows R in 
blue. We will use the renormalization scheme to classify the periodic points 
of*. 
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Lemma 7.2 Let 0\ = 0~ fe (0i) and 0\ = 0" fc (O 2 ). Then 

• There is a self-similar packing of regular octagons in S , each of which 
is translation equivalent to 0\ (or 0\) for some k. (See Figure 7-4-) 

• Every periodic point of $ belongs to one of the octagon interiors, and 
its period is the same as the tile-period of the octagon containing it. 

• Every octagon in the packing has tile period 3 fc for some k — 0, 1, 2, 
and such an octagon either lies in the orbit of 0\ or the orbit of 0\. 
These are the blue octagons in Figure 7.4- 

Proof: We see by inspection that <j) 3 (S') C S'. Indeed, the left hand side of 
Figure 7.5 shows S in blue, the middle picture shows ^(S) in blue, and the 
right hand side shows $ 2 (S') in blue. 




Figure 7.5: Image of the central portion 



Lemma 7.3 Every ^ -periodic point in R has period 3 fc for some k = 0, 1, 2... 

Proof: A short calculation establishes the following two facts. 

0- 1 o$o6| 5 = ^ 3 ; R-0 1 -0 2 = SU^(S)U^ 2 (S). (35) 

The same analysis as in the previous chapter produces the packing of oc- 
tagons, together with the claims about their tile periods. The only new 
feature of the situation here is that the map ^ is a piecewise translation, 
whereas the map from the previous chapter had a rotational component. 
Since \1/ is a piecewise translation, each point in a octagonal tile has exactly 
the same period as the tile itself. 4k 
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7.5 Proof of Lemma 15.11 



Thanks to the analysis done in £J5] it suffices to prove that every periodic 
point of R\ has $-period 3 k for some k = 0,1,2, .... We have already worked 
out the structure of the periodic points of the \l/-system, and we will use this 
to get the desired information about the periodic points of the $-system. 

There is a 2 to 1 map F : Ri — y R. which, by construction, carries 
periodic points of $ to periodic points of Moreover, this map conjugates 
the action of $ on R l to the action of \1/ on R. From this, we see that every 
^-periodic orbit has period either 3 fc or 2 x 3 fe . 

A parity check rules out the second possibility. Given a \I/-periodic orbit 
Po,Pi, ■■■jPn-i, we let crj = if p lies in a blue region of R and other wise 
<jj = 1. We just need to check that Ylu=o a i ls even - We check that the 
parity of the sum is preserved under the renormalization operation. Since 
the sum obviously works out for the large octagons - i.e., when iV = 1, it 
works out in general. 

The parity check is a finite calculation, as we now explain. Suppose we 
start with a point z G R. Let w = 9 _1 (i?). Let Wj = ^(wo). We just have 
to prove that 

cr = cr + <j 1 + (T 2 mod 2. (36) 

Given the nature of we only need to check this equation for 1 point in 
each of the 22 regions. Here are the results for the first 11 regions 

0:0 111 

1:0 000 

2:1 001 

3:0 000 

4:1 100 

5:0 101 

6:1 111 

7: 1 111 

8:1 111 

9:1 111 

10 : 10 1 

The pattern repeats exactly for the second 11 regions. 
This completes the proof of Lemma 15.11 
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8 Generating the Arithmetic Graph 



8.1 The Basic Approach 

Let R\ be our basic region, partitioned as on the left hand side of Figure 7.1. 
To each region r 6 R\ we choose and pgr and define 

<fc(r)=7r fc (S(p)). (37) 

Here € -Z 8 is the vector from Equation |H Since $ is constant on the 
region r, we see that our definition is independent of the choice of p G r. If 
p is part of some orbit, then gk(p) € R 2 is the side of the edge of TCk(T) that 
corresponds to (p, 

Given an initial point po £ -Ri we consider the forward orbit po,pi,p2, ■■■ 
We have a sequence of regions r ,ri,r 2 ,... and a corresponding sequence 
9k(ro) , gk(ri) , g k (r 2 ) , ... These vectors are the successive sides of the projection 
TTfc o r(p). 

8.2 The Compressed Approach 

Given the nice renormalization scheme we discussed in the previous chapter, 
we prefer to work with the compressed system \1/ : R — > R, as defined in §7.31 
The region R is colored as in Figure 7.3. 

Recall that R is the top half of R\. Let R' be the bottom half. Let 
Po,pi,P2, ■■■ be a f orbit in R. Let go, <?i, <?2 be a $ orbit in Ri such that 
Po = go- Let Oi = if pi lies in a blue region of R on the left hand side of 
Figure 7.3, and (Tj = 1 otherwise. 

We define the accumulated parity 

t=Q 

That is, we count the number of times that the orbit lands in a red region 
on the left hand side of Figure 7.3. The sequence {s n } is a binary sequence 
that keeps track of whether p n = q n or not. 

The accumulated parity is the additional information that lets us "lift" 
an \1/ orbit on R to a $ orbit on R±. 

Lemma 8.1 s n = if and only if p n = q n . 
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Proof: This follows from two facts. First, Oi = iff and qi lie in the 

same half of Ri, Second, $ : R\ — > R\ commutes with the map interchanging 
R and R' . 6 

To each region r G R, we associate 4 vectors according to the following 
rule. 

9k(r,0) = g k (r); g k (r, 1) = g k (r'); r' = r - (0, 1) e R'. (39) 

Given the orbit Po,Pi,P2 G -R we have the corresponding regions r ,ri,r 2 , ... 
and the corresponding accumulated parities so, si, S2, By construction, 
^fc(r n , s n ) is the nth edge of 7r fc r(p )- 

For reference, we list the vector assignments. It is convenient to let k 
stand for the (region, parity) pair (k, 0) and let k + 22 stand for the (re- 
gion,parity) pair (A;, 1). The indices 0,11,22,33 correspond to the biggest 
octagons and never arise, so we leave them off [HI We consider the gi vectors 
first. Letting T(x, y) = (—y, —x), we have the symmetries 

g 2 (k + 11) =T(g 2 (k)); g 2 (k + 22) = —g 2 (h)\ fc = 0,...,10. (40) 

So, the following assignment determines all 44 vectors. 



1. 


ffa(i) = 


(6,2). 


2. 


02(2) = 


(2,-2) 


3. 


02(3) = 


(2,-2) 


4. 


02(4) = 


(-4,4) 


5. 


02(5) = 


(-6,2) 


6. 


2 (6) = 


(0, -4) 


7- <fe(7) = 


(2,6). 


8. 


2 (8) = 


(4,4). 


9. 


02(9) = 


(2,2). 



10. 2 (1O) = (-2, 6). 

3 For the record, g 2 (0) = (0,8) and g 3 (0) = (0,0). 
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In the case of g 3 we have the symmetries 

g 3 (k + 22) = g 3 (k); g 3 (k + 11) = -g 3 (k) +a k (0, 4); k — 0, 10. (41) 

Here <7 fc G {0, 1} is the parity of the kth region. We have 

a k = l =}► fc e {2,4,6,7,8,9}. (42) 

In light of these symmetries, the following list of 10 vector assignments de- 
termines the whole thing. 

1. g 3 (l) = (0,-1, -2,-1). 

2. <7 3 (2) = (2, 1,0,-1). 

3. «? 3 (3) = (2, 1,0,-1). 

4. (? 3 (4) = (-2,-2,2,0). 

5. <7 3 (5) = (-2, -1,0,-1). 

6. ^ 3 (6) = (0,0,0,-2). 

7. 5s(7) = (-2,-1,4,1). 

8. <fe(8) = (-2,-2,2,0). 

9. 3 (9) = (-2,-1,0,-1). 
10. </ 3 (10) = (0,-1,2,1). 

Our notation is as follows. 

(a, 6, c, d) = (a + bV2,c + dV2). (43) 

In the next section we will explain how these vector assignments combine 
with a combinatorial substitution rule to generate the graphs. 
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8.3 The Substitution Approach 



Suppose that we have some orbit Pq,Pi,P2, ■■■■ We let 

Po = ®~\Po) 

and consider the renormalized orbit p' ,p'i,P2, • Let Sj = <To + ... + <Tj-i be 
the accumulated parity of the original orbit. Likewise define s'j. Let rj be 
the region of R containing pj, and likewise define r'j. 

In §7.51 we explained the sense in which the renormalization scheme re- 
spects the accumulated parities. This fact implies that 



Given equation [44], the pair (sj,rj) determines the triple 

( r 3i' S 3j) ( r 3j+l' S 3j+l) ( r 3j+2' S 3j+2)- 

Recall that we have chosen a particular listing for these 44 possible (region, 
parity) pairs. What we are saying is that the numerical code for the original 
orbit determines the numerical code for the renormalized orbit, by way of a 
substitution scheme in which one number is replaced by 3 numbers. Here is 
one quarter of the substitution rule. 






— >■ 


9 


25 


39 


1 


-> 


10 


14 


5 


2 


— >■ 


10 


12 


17 


3 


->■ 


10 


12 


16 


4 


-> 


9 


23 


27 


5 


-)■ 


8 


36 


28 


6 


-> 


8 


35 


17 


7 


-> 


9 


24 


6 


8 


-> 


9 


26 


6 


9 


—> 


9 


26 


7 


10 


-> 


9 


25 


40 



The rest of the pattern can be deduced from symmetry and from what we 
have already written. We point out two involutions on the set {0, 43}. The 
first of these has the action n — > n ± 22. The second one has the action n — > 
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n ± 11, where the + option is taken if n e {0, 10} or n e {0, 10} + 22. 
The substitution rule is invariant under these two involutions. For instance, 

4 -> 9 23 27 

15 -»■ 20 34 38 

26 4 31 1 5 ' 

37 -»■ 42 12 16 

Now we can explain a nice combinatorial way to generate the arithmetic 
graph projections. Referring to Lemma I7.2[ the two orbits 0\ and 0\ both 
have period 3 fc . Here 0\ and 0\ respective are the top and bottom biggest 
octagons in Figure 7.4. We will concentrate on the orbits 0\. The other half 
of the orbits have a very similar treatment. 

We start with 0, and iterate the substitution rule k times. Then we 
replace each number by the relevant vector. The resulting list of vectors 
is the set of edges of the projection of the arithmetic graph. We call the 
resulting path G 2 {k) or G 3 (k), depending on whether we use the G 2 vectors 
or the G3 vectors. 

Here are two examples. For G 2 (2) we make the substitution — ?- 9, 25, 39, 
and then plug in the G 2 vectors: 

• 94 (2,2). 

• 25 4 (-2,2). 

• 39 ^ (0,-4). 

We then get the vertices of G 2 (l) by accumulating the vectors: 

(0,0), (2, 2), (0,4), (0,0). 

So, G 2 (l) is a closed path - a triangle. 
To generate G 2 (2) we go one more step. 

0^ 9, 25, 39 4 9, 26, 7, 32, 34, 38, 41, 2, 28. 

Then we make all the substitutions and accumulate the vectors. We discover 
in this case that G 2 (2) is not closed. We check that G 2 (k) is closed if k = 
1,3,5,7 and open if k = 2,4,6,8. In §9.31 we will prove that this pattern 
persists for all k. Our proof there, in particular, implies Statement 1 of the 
Main Theorem. 
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9 Taking a Better Limit 



9.1 Subsequential Limits 

Note that G 2 (k) and G 3 {k) have natural parametrizations, coming from the 
ordering on the edges. Define 

T 2 (k) = (V3)- k G 2 (k)- F 3 (k) = (1 + V2)- k G 3 (k). (45) 

We notice that T 2 (k) and T 2 (k + (3) closely follow each other as parametrized 
curves when (3 = 4. This property does not hold for smaller (3. When (3 = 2, 
the two curves T 2 (k) and T 2 {k + (3) are close in the Hausdorff metric, but 
not close as parametrized curves. We also notice that T 3 (k) and T 3 (k + (3) 
closely follow each other closely when (3 = 2 but not when (3 = 1. 

Note that (3 = 4 is the lowest number that works well for both curve 
families. This accounts for the dependence on the mod 4 congruence that 
appears in the Main Theorem, which involves the subsequential limits 

limr i (i + 4A;); * = 2,3; j = 1,3. (46) 

The 4th power of our substitution rule is our basic object of study. With 
this rule, each number is replaced by 81 numbers. In geometric terms, each 
G 2 vector g 2 (n) is replaced by a polygonal path g' 2 (n) of length 81. This path 
consists of 81 other G 2 vectors, scaled down by \^3. Likewise, each G 3 vector 
g 3 {n) is replaced by a polygonal path g' 3 (n) of length 81. 

Let Rg 2 (n) denote the vector that points from the tail of g' 2 {n) to the head 
of g' 2 {n). Likewise define Rg 3 {n). We notice that gi(n) and Rgiin) are nearly 
the same vector for each choice % = 2, 3 and each n = 0, ...,43. We think of 
R as a renormalization operator, which replaces one vector by another one 
that is very close to it. It makes sense to define R k gi{n) as follows: 

• We iterate the substitution rule Ak times. 

• We replace each number by the appropriate ^ vector. 

• We rescale by (si)~ k , where s 2 = \/3 and s 3 = 1 + \[2. 

The key to analyzing our limits is to looked at the fixed point of the renor- 
malization operator. 
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9.2 The Fixed Point of Renormalization 



For motivational purposes (only) we mention the Perron-Frobenius Theorem. 

Theorem 9.1 (Perron-Frobenius) Let T : R n — > R n be a linear trans- 
formation represented by a matrix with all positive entries. Then T has ( up 
to scale) a unique unit positive eigenvector u. Moreover, if v is any positive 
vector v, then we have v n — > u, where v n = T n {v) /\\T n {v)\\. 

It turns out that gi(n) and Rg,i{n) are not in general the same vector. 
The Perron-Frobenius Theorem suggests to us that we consider the alternate 
assignment 

n ->• Xi(n) := lim R s gi(n). (47) 

s— >oo 

Doing things this way is analogous to finding the Perron-Frobenius eigen- 
vector of a transformation T by taking the limit of iterates of T applied to 
some conveniently chosen initial vector v. In case the limit u is expected to 
have nice properties - e.g. algebraic coordinates - one might hope to guess u 
by looking at the decimal expansions of a fairly high iterate. Once we have 
made a guess, it is trivial to verify that it works. This is what we did. 

Our new assignment for G 2 has the same symmeties as in Equation HOJ 
so the 10 vectors we list suffices to describe the whole thing. 



1. 


A 2 (l) = 


* (6,2). 


2. 


A 2 (2) = 


(3,-3). 


3. 


A 2 (3) = 


* (2,-2). 


4. 


A 2 (4) = 


(-3,3). 


5. 


A 2 (5) = 


* (-6,2). 


6. 


A 2 (6) = 


(-3,-1) 


7. 


A 2 (7) = 


(3,5). 


8. 


A 2 (8) = 


(5,3). 


9. 


A 2 (9) = 


(3,1). 


10. 


A 2 (10) : 


=* (-2,6) 
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This assignment has the same symmetries as the original one, and these 
symmetries let one deduce the remaining vectors in the assignment. The 
starred equalities indicate that no change has been made. Curiously, we 
have 

new = old + e(l,-l); eG {-1,0,1}. (48) 

In the case of G3, the new assignment has nicer symmetries than the 
original namely 

A 3 (A0 = -A 3 (A; + 11) = A3(M). (49) 

That is, the factor of is replaced by in Equation [JTJ In light of these 
symmetries, our 10 vectors below determine the whole thing. 



1. 


A 3 (l) = 


* (0,-1,-2,- 


2. 


A 3 (2) = 


(2,1,-2,-1). 


3. 


As(3) = 


* (2,1,0,-1). 


4. 


A 3 (4) = 


(-2,-2,0,0). 


5. 


A 3 (5) = 


* (-2,-1,0,- 


6. 


A 3 (6) = 


(0,0,-2,-2). 


7. 


A 3 (7) = 


(-2,-1,2,1). 


8. 


A 3 (8) = 


(-2,-2,0,0). 


9. 


A 3 (9) = 


(-2,-1,-2,- 


10. 


A 3 (10) : 


=* (0,-1,2,1). 



Again, the starred equalities indicate the ones where there is no change. This 
time we have 

new = old + e(0,2); eG {-1,0,1}. (50) 

A direct calculation shows that 

R\i(n) = Ai(n); i = 2,3; ra = l,...,43. (51) 

As usual, we are omitting n = 0, 11, 22, 33. 
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Let Aj(fc) denote the version of T^k) produced by the new vector assign- 
ments. Our new assignments are nicer, but we need to prove that the A 
paths have the same limits as the T paths. Once we know this, we can throw 
out the arithmetic graphs, and work with the much nicer A paths. The rest 
of the chapter is devoted to analyzing and equating the relevant limits. 

Figure 9.2 shows r 3 (3) and r 3 (7). Notice that these curves follow each 
other pretty well, but that a few of the vertices of T 3 (7) do not lie on T 3 (3). 
Figure 9.3 shows the improved A 3 (3) and A 3 (7). These curves are better 
related to each other. Note also that T 3 (7) and A 3 (7) are really just about 
the same curve. 




Figure 9.2: T 3 (3) and T 3 (7). 
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Figure 9.3: A 3 (3) and A 3 (7). 



9.3 Proof of Statement 1 

Statement 1 of the Main Theorem, which only deals with the odd case, follows 
from the more complete result we now prove. 

Lemma 9.2 The following is true. 

• Gs(k) is a closed polygon for all k = 1, 2, 3, ... 

• G 2 (k) is a closed polygon if and only if k is odd. 
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Proof: Let be the unsealed version of A*.. The main idea in our proof is 
showing that Gk and have the same endpoints. Assuming this result, it 
suffices to prove our result for the L curves. But the A curves are just scaled 
version of the L curves. So, it suffices to prove our result for the A curves. 
Given the endpoint-preserving property of our ideal substitution rule - see 
Equation ED - we see that Aj(fc) is closed if and only if Ai(k + 4) is closed. 
We check the conclusion of this lemma for the first few cases, and then apply 
the induction step to take care of the remaining cases. 

Now we come to the interesting part of the proof, showing that each G 
curve has the same endpoints as the corresponding L curve. We treat G2 
and L 2 first. All paths start at the origin, so we just have to see that the 
other endpoints match. That is, we need to prove that 

AT-l 

S(rn) = 0; S(m) = g 2 (rii) - A 2 (n i ). (52) 

Here {rii} is the numerical code for the relevant pair of paths, and N is the 
total combinatorial length. 

As we remarked above - see Equation|3H]~ we always have S(n) = e(l, — 1) 
for eG {-1,0,1}. We check by hand the following 44 equations. 

6{n) = 6{m 1 ) + 6{m 2 ) + S{m 3 ). (53) 

Here (m 1 ,m 2 ,m 3 ) is the substitution for n. Equation |5"31 implies that the 
sum in Equation [52] is invariant under substitution. We check that the sum 
is in the first 4 cases, and then we see by induction that it is always 0. 

The proof for G3 and L 3 works in the same way. This time, we have 
5{n) = e(0, 2) for e G {-1, 0, 1}. * 



9.4 Proof of Statement 2 

Let G and L be the unsealed versions of V and A respectively, as in the proof 
of Lemma 19.21 Say that a distinguished strand of G is one that comes from 
repeated substitution applied to a single edge. Likewise define distinguished 
strands of L. Such strands consist of 3 fc segments. Any distinguished strand 
a determines a vector V(a) that points from the starting point to the ending 
point. 
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Lemma 9.3 Let g and I be corresponding distinguished strands of G(n) and 
L(n). Then V(g) and V(l) differ by at most 2 units. 

Proof: Let g' and I' be the distinguished strands of the paths G(n — 1) and 
L(n — 1) which give rise to g and /, respectively, upon substitution. The same 
argument as in Lemma 19.21 shows that 

V(l')-V(g')=V(l)-V(g). 

By induction, we are reduced to the case of single edges. Now we apply 
Equations [48] and [50] 4k 

Corollary 9.4 The distance between any corresponding vertices of Gi(n) 
and Li{n) is at most 41og 3 (n). 

Proof: Both Gi{n) and Li(n) start at (0,0). Suppose we want to consider 
the kth vertices. Let 7 be the arc of Gi connecting (0, 0) to the fcth vertex. 
Likewise define A. 

Using the base 3 expansion of A; as a guide, we can decompose 7 and A 
into at most 21og 3 (n) distinguished arcs. Our lemma now follows from the 
previous result and the triangle inequality. 4ft 

Lemma 9.5 Let i be 2 or 3. The sequences {Aj(l + 4fc)} and {Aj(3 + 4fc)} ; 
considered as a sequence of parametrized curves, converge in the L^ norm 
as k 00. 

Proof: We will consider the first of these sequences. The second one has the 
same proof. The edges of Aj(4/c + 5) are obtained by replacing each edge of 
A(4/c + 1) by a length 81 polygonal path that has the same endpoint. Define 

D k = sup dist(Ai(4£; + l;s),Ai(4£; + 5;s).Y (54) 
se[o,i] ^ ' 

Each edge of A(4A; + 5) is 

1/9 = (V3)- 4 
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times as long as an edge of the same type in A(4k + 1). By similarity, we 
have 

D k+1 = -±=D k , (55) 

once k is large enough that Aj(l + 4k) contains all the edge types. It suffices 
to take k = 1, in fact. 

Equation [55] implies that the distance between two successive curves in 
our sequence decreases exponentially. Hence, our curves form a Cauchy se- 
quence in the space Map(S' 1 , -R 2 ), equipped with the topology. Hence, 
our sequence converges. 4fc 

The following result implies Statement 2 of the Main Theorem. 

Lemma 9.6 The sequences {Tj(l + 4k)} and {Tj(3 + 4k)}, considered as a 
sequence of parametrized curves, converge in the norm as k — > oo. The 
limits in this case coincide with the limits in Lemma \9. 51 

Proof: We prove this first for i — 2. The i — 3 case is essentially the same. 
By Corollary l9.4l the distance between corresponding points of G2(l+4k) and 
L 2 (l + 4k) is at most 41og(8A;). But then, by scaling, the distance between 
corresponding points of T 2 (l+4k) and A 2 (l+4/c) is at most 4 \og 3 (8k)(V3)~ k , 
a quantity that tends to exponentially fast. 4 
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10 The End of the Proof 



10.1 Scaling Constants 

In this chapter, we prove Statements 3, 4, and 5 of the Main Theorem. 

Let I 2 and J 3 be suitably scaled versions of the sets discussed in §HJ For 
i — 2,3 and j — 1,3 let Ay be the limit of the sequence {Aj(j + 4fc)}. To 
finish the proof of the Main Theorem, we need to show that 

A 21 = A 23 = h; A 31 = A 33 = I 3 . (56) 

Our first order of business is to explain how to scale I 2 and I 3 . 




Figure 10.1: J 2 seed and A 2 (l). 



We will deal with J 2 first. Figure 10.1 shows the initial seed for J 2 , scaled 
so that it is well situated with respect to A 2 (l). By construction, A 2 (l) is 
an isosceles triangle with vertices (0,0) and (3, 1) and (1,3). Given the way 
that the J 2 seed sits with respect to this triangle, we deduce that the center 
of I2 is (3/2,3/2) and the left corner is (—3/2,3/2). This gives us the copy 
of I2 mentioned in the Main Theorem. 
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Now we deal with J 3 . Figure 10.2 shows A 2 (l) and the union i 3 (l) of 
triangles obtained by applying the substitution operator to the I 3 seed. 




Figure 10.1: J 3 (l) and A 3 (l). 



A 3 (l) is an isosceles triangle, with vertices 

(0,0); (_2 - >/2, -2 - >/2) (0,-2-2^2). (57) 

The point (0, 0) is the apex of the triangle, the top vertex. 

This is enough to determine the placement of J 3 . Here is the nicest 
way we can think of describing the scaling. Let s = 1 + y/2. This is the 
scaling factor that came up in the previous chapter. The I 3 seed consists 
in 8 isosceles triangles. One of the isosceles triangles has vertices (0,0), 
(1/2, — s/2) and (— s/2, —1/2). This gives us the version of J 3 mentioned in 
the Main Theorem. Note that all the vertices of J 3 , as we have scaled it, lie 
in Q[V2\. 
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10.2 Proof of Statement 3 



Let I2 and ^3 be scaled as in the previous section. By Equation [5H we see that 
Ay contains the vertices of all the polygons of which it is a limit. Moreover, 
by scaling, the vertices of the approximating polygons become dense in 
as the number of sides tends to 00. We passed to the limit of renormalization 
precisely to arrange these two properties. We want to understand the sets in 
Equation [56] just from the placement various finite collections of vertices. 

Given the properties enjoyed by the vertices, we can establish Equation 
[561 simply by showing, for % = 1, 2 and for all odd j that the vertices of the 
polygons approximating Ay are contained in Jj and become dense in as 
j -> 00. 

We will deal with % = 3 first. Figure 10.3 shows a piece of the 13(3), 
together with a piece of A3 (3). Notice that all the vertices of A3 (3) in sight 
coincide with right-angled vertices of is(3). Notice also that there seems to 
be (most of) a smaller copy of Figure 10.1 sitting near the top of Figure 10.3. 




Figure 10.3: J 3 (l) and A 3 (l). 
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Say that a special vertex of J3 is a vertex of a right-angled triangle of 
some .^(n). We check by direct (computer) calculation that every vertex of 
A^(j) is a special vertex of -Z3 for j = 1, 3, 5, 7, 9. This is an easy calculation 
that takes place entirely inside Q[y2). 

We can't make an infinite number of explicit calculations to deal with 
each successive polygon, so we describe a method, based on self-similarity, 
that proves everything we need from a finite number of calculations. For 
ease of exposition, we will work entirely with the sequence A 3 (l + 4k). The 
sequence A 3 (3 + 4/c) has a similar treatment. It is convenient to set A 3 = A 31 . 
So, A 3 is the limit of the sequence we are considering. 

The polygon A 3 (l) is just a triangle, and too small to be of any use 
to us (beyond what we have already done with it above.) So, we consider 
the next polygon in the sequence. Each edge of A 3 (5) is a scaled copy of 
some A 3 (n). We call n the type of the edge. The type belongs to the set 
S = {1, ...,43} — {11,22,33}. Inspecting the numerical sequence associated 
to A 3 (5) - see §8.31 - we check that A 3 (5) has every possible type edge. We 
check the same thing for A 3 (9). 

Suppose that A C A 3 (5) is an edge of type n and B C A 3 (9) is an edge 
of type n. The pair (A, B) determines a homothety 

H = H(A,B) : A^ B. (58) 

Given the way A 3 is the limit of our sequence, we have 

H:A 3 (A)^A 3 (B). (59) 

Here A3 (A) is the portion of A 3 between the endpoints of A, and likewise for 
A 3 (B). So, if is a homothety that carries a small portion of A 3 to an even 
smaller portion. We call H a special homothety associated to the pair (A, B). 
Every such pair gives rise to a special homothety. We call H good if 

H:h{A)^h{B). (60) 

Here 13(A) is the portion of I3 bounded by the endpoints of A. Likewise we 
define 13(B). This definition makes sense because we have shown that the 
endpoints of A and B both lie in J3. Just to be clear, we mean that H sets 
up a bijection between 13(A) and 13(B) that maps endpoints to endpoints. 
Below, we will prove the following result. 

Lemma 10.1 (Good Homothety) Every special homothety is good. 
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The Good Homothety Lemma involves a finite number of special homotheties. 
We will show in the next section how one verifies that a special homothety 
is good using a finite calculation. Indeed, one can practically see this by 
inspection using our program OctoMap 2. 

Before we prove the Good Homothety Lemma, let us explain what it does 
for us. Suppose we want to show that the vertices of A 3 (13) lie in J 3 . Let v 
be a vertex of A3 (13). If v is also a vertex of A3 (9), we are already finished. 
Otherwise, there is some edge B of A3 (9) such that v is one of the vertices 
associated to the replacement of B. That is, v lies on the polygonal path of 
A 3 (13) connecting the endpoints of B. Choose and edge A of A 3 (5) that has 
the same type as B. 

The special homothety H = H(A, B) carries A to B. Moreover, H carries 
A 3 (9; A) to A 3 (13; B). Here A 3 (9; A) denotes the portion of A 3 (9) subtended 
by A. At the same time as this, H carries 13(A) to 13(B). But H^ Y (v) is a 
vertex of A3 (9; A). Hence H~ x (v) 6 h(A). Hence v G h(B). This proves 
that all vertices of A 3 (13) lie in I3. 

Continuing in this way, we see that all vertices of A 3 (l + 4/c) lie in 7 3 for all 
positive k. To take care of density, we let c4 denote the maximum distance 
from a point of I3 to a vertex of A 3 (l + 4k). Since our special homotheties 
scale distances by a factor of \/2 — 1, we see that 

4+i = (\/2-i)4- (61) 

In particular dk — > as k — > 00. This proves that the vertices of A 3 (l + 4k) 
becomes dense in 1 3 as k — > 00. 

In short, the Good Homothety Lemma implies Statement 3 of the Main 
Theorem. 

10.3 Proof of the Good Homothety Lemma 

In §4.21 we described the basic substitution rule for the snowflake, in which 
each isosceles triangle is replaced by 5 smaller ones. We started with a 
seed, consisting 8 isosceles triangles arranged in a cyclic pattern, and then 
iteratively replaced each triangle by the smaller ones. 

Here we describe a more conservative approach to the construction. Let 
L\ and L 2 be finite cyclically ordered lists of isosceles triangles. We write 
Li — > L 2 if one obtains L 2 from L\ by replacing one of the triangles of L 2 by 
the 5 smaller ones from the substition rule and adjusting the cyclic ordering 
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in the obvious way. We call a list L of triangles admissible if we have a finite 
chain ^(O) = L — > ... — > L n . The specific lists h{n), produced in §4.21 are 
examples of admissible lists. 

We say that a special vertex of I3 is a right-angled vertex of some Is(n). 
An equivalent definition is that a special vertex of ^3 is a right-angled vertex 
of some admissible list. Let {v k } be a finite list of special vertex. We say that 
an admissible list L is compatible with {v k} if each point v k is a vertex of some 
triangle in L. We say that L is the minimal pattern for {vk} if L is compatible 
with {vk} and there is no V such that V — > L and L' is compatible with {vk}. 

Remark: For our proof of the Good Homotopy Lemma, we do not need 
to know that the minimal pattern associated to a finite list of vertices is 
unique, but we emphasize that the minimial pattern is unique. The unique- 
ness comes from the tree-like nature of the subdivision rule. We simply start 
with the seed and replace triangles only when necessary. After a finite num- 
ber of steps we arrive at the minimial pattern. 

Let A be a line segment whose endpoints are special vertices of I3. We 
let Js(A) be the minimal pattern associated to the endpoints of A. We don't 
care so much about the whole pattern, but only the portion related to A. We 
let Jz(A) denote the portion of J^(A) that lies between the two endpoints of 
A. Figure 10.4 shows the minimal pattern J${A) when A is a certain type-1 
edge of A 3 (5). 



Figure 10.4: Jz(A) for A a type 1 edge of A 3 (5). 
The figure we have shown is completely representative of the local picture 




78 



for the type 1 edges of A 3 (5). The associated minimal patterns are all the 
same. This is not really so surprising, because these vectors all are translates 
of each other. We also check that the local picture looks the same for all type 
1 edges of A 3 (9). That is, the picture looks the same up to a scale factor of 



So, suppose that A is a type 1 edge of A 3 (5) and B is a type 1 edge 
of A 3 (9). Then H(A,B) carries Jz{A) to J%{B) because the patterns are 
identical up to scale. But J%{A) and J%{B) respectively determine / 3 (^4) and 
hiB). Hence H carries 13(A) to h(B), as desired. This proves the Good 
Homotopy Lemma for edges of type 1. The key idea is that the minimal 
pattern J 3 (A) only depends on the type of A, up to scale. 



Figure 10.5: The minimal patterns associated to types 2,3,4. 

Figure 10.5 shows the minimal patterns associated to the types 2, 3, and 
4. All the edge types except those of the form 6 + Ilk and 9 + Ilk have 
the same property we have already discussed. The local picture is the same, 
independent of the edge. By inspection, we see that the Good Homotopy 
Lemma holds for for all edge types except the 8 types we have mentioned. 

For each of the remaining edge types, there are two local pictures, up to 
rotation, these two local pictures are the same for each of the edge types. 
Figure 10.6 shows the two pictures associated to the edges of type 6 + Ilk. 
The edges are vertical. 
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Figure 10.6: The minimal patterns associated to the edges of type 6 + Ilk. 

Let A and A' be the two type 6 edges shown in Figure 10.6. The pattern 
J$(A) consists of 4 triangles T1UT2UT3UT4 and J$(A') consists of 3 triangles 
T[ U T' 2 U T3. We list these triangles top to bottom. There is a translation r 
such that t(A) = A'. To finish our proof for edges of type 6 + Ilk, we just 
have to see that r carries 13(A) to I 3 (A'). Note that r carries T- to Tj for 
j = 1,2. So, things work out as desired for the "top portions" of the sets of 
interest to us. 

Now we consider the bottom portions. Let S be the line segment on the 
left hand side of Figure 10.6 that connects the bottom vertex of A to the 
bottom vertex of T 2 . Note that S is made from the hypotenuse of T 3 and a 
short side of T 4 . Let S' be the edge of T 3 that connects the bottom of A' to 
a vertex of T' 2 . We see by direct calculation that r(S) = S'. The point is 
that S and S' have the same slope and length. We just have to prove that r 
carries h(S) to h(S'). But we have exactly the situation discussed in §4.31 
and the desired result here follows from Lemma 14.11 

it remains to consider The edges of type 9 + Ilk. It turns out that 
the minimal patterns associated to these edges look exactly like the ones 
associated to the edges of type 6 + Ilk, except that the picture has been 
rotated 45 degrees. The same argument as in the previous case works here. 

This completes the proof of the Good Homothety Theorem. Again, the 
Good Homothety Theorem (and our previous work) implies Statement 3 of 
the Main Theorem 
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10.4 Proof of Statement 4 



In considering J 2 we make all the same definitions that we made for I 3 , 
except that we redefine what we mean by a special point. In the case of I3, 
the special points are the right-angled vertices of the triangles. In the case 
of I 2 , we call a point special if it has one of two properties: 

• It is the center point of a parallelogram of some admissible list. 

• It is the center of the long edge of a trapezoid of some admissible list. 

We check directly that every vertex of I 2 (j) is a special point, for the initial 
values j = 1, 3, 5, 7, 9. This check entirely involves calculations with rational 
points. We want to extend this result for all odd j. As in the A3 case, we 
ignore the sequence A(4k + 3) and concentrate on the sequence A(4k + 1). 
The other sequence has essentially the same treatment. 

To finish the proof, we just have to prove the analogue of the Good 
Homothety Lemma for I 2 . Given the combinatorially identical subdivision 
rules that produce I2 and I3, it makes sense to say that a minimal pattern 
associated to I2 is isomorphic to a minimal pattern associated to J 3 : The 
canonical bijection between the /3-shapes and the /3-shapes carries the one 
pattern to the other. For j = 1, 3, 5, 7, 9 we make the following observation: 
The minimial pattern associated to some edge of A2(j) is isomorphic to the 
minimal pattern associated to the corresponding edge of A^(j). There is a 
perfect correspondence between the two kinds of patterns. 

Given the combinatorial isomorphism between the /3-patterns and the 
7 2 -patterns, the proof of the Good Homothety Lemma goes through, word 
for word for all the edges except those of type 6 + Ilk and 9 + Ilk. In these 
cases, there are two combinatorial types of minimial pattern, and we must 
make a separarate analysis. Figure 10.7 shows the relevant portions of the 
minimal patterns associated to edges of type 6 + Ilk. 
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Figure 10.7: The minimal patterns associated to the edges of type 6 + Ilk. 

On the left hand side we have a chain of 3 quadrilaterals and on the right 
hand side we have a chain of 4 quadrilaterals. The first two pieces agree on 
both sides. When we chop two pieces off, and superimpose what remains of 
the left and right sides, we get exactly the left hand side of Figure 4.7, the 
figure we used to illustrate the hidden symmetry of the carpet. The rest of 
the proof is just like what we did for the A 3 case, except that we use Lemma 
I4.6t n place of Lemma 14.11 

The edges of type 9 + Ilk have the same treatment. This completes the 
proof of the Good Homotopy Lemma in the I2 case, and thereby completes 
the proof of Statement 4 of the Main Theorem. 

10.5 Proof of Statement 5 

Statement 5 takes a bit of unwrapping. We start with an arbitrary sequence 
{x n } of good points in the half-strip Ej, and we let A n be the arithmetic 
graph corresponding to x n . We want to understand the map 



Here S^n and Ss >n are the dilations discussed in the Main Theorem. 

In light of the analysis done in the previous chapter, we can assume that 



4>n ■ S; 



3,^3 (Ai) -> S 2 ,n^2{A n ) 



(62) 



3, 



'k,n^k{A n ) = r k (n). 



(63) 
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Here n is odd, and T^n) is as in §9.11 For ease of exposition, we now make 
the blanket assumption that n = 1 mod 4. The case when n = 3 mod 4 has 
the same treatment. 

The two paths r 2 (n) and r 3 (ri) have the same number of sides, and so 
there is a canonical map 

n : r 3 (n) -> r 2 (n). (64) 

There is a canonical map from ^(n) to Afc(n), and this gives us a map 

V> n :A 3 (n)->A 2 (n). (65) 

Our strategy is to first analyze the sequence {ip n } of maps and then use 
Lemma [9751 to get the desired information about {<p n }- 

There is a sense in which ip n and <fi are compatible, we we now explain. 
The map <p is the limit of a sequence of canonical bijections between the 
admissible J 3 -patterns of triangles and the admissible J 2 -patterns of quadri- 
laterals. We call this bijection the shape bijection. We denote it by $. Note 
that $ is canonical only up to composition with isometries of / 2 and J3. 
This accounts for the isometry F that appears in Statement 5 of the Main 
Theorem. 

Given an edge e 3 of A 3 (n), there are two triangles 7\ and T 2 in the minimal 
pattern A 3 (e 3 ) whose right-angled vertices are the endpoints of e 3 . Let 

e 2 = ^n(e 3 ) 

be the edge of A 2 corresponding to e 3 . We say that the pair (e 3 ,e 2 ) is 
compatible if the endpoints of e 2 are special points of the two shapes $(Ti) 
and $(T 2 ). In case $(Xi) is a parallelogram, the relevant endpoint of e 2 
should be the center point of this parallelogram. In case <&(Ti) is a trapezoid, 
the relevant endpoint of e 2 should be the midpoint of the long side of this 
trapezoid. 

We check that all relevant edge-pairs are compatible for n — 1, 3, 5, 7, 9. 
But our proofs of the Good Homotopy Lemma, in the two cases, involved 
combinatorially identical sets in each case. That is, whatever we did in the 
Is case carries over to the 7 2 case via the shape bijection. From this we see, 
by induction and symmetry, that all relevant edge-pairs are compatible for 
all odd positive integers n. 

Given this notation of compatibility, we have 

lim ip n = F o 0, (66) 
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for some isometry F. The convergence takes place in the sense that the graph 
of i/) n converges, to the graph of in the Hausdorff topology. In particular, 
the left hand side of Equation [66] does have a limit. But Lemma [9T61 now says 
that 

lim 4> n = lim ijj n . (67) 

Putting the two equations together gives us Statement 5 of the Main Theo- 
rem. 

This completes our proof of the Main Theorem. 
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